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Abstract. We investigate the formation of the ten double-lined double white dwarfs that have been observed so far. A detailed 
stellar evolution code is used to calculate grids of single-star and binary models and we use these to reconstruct possible evolu- 
ti onary scenarios. We a pply various criteria to select the acceptable solutions from these scenarios. We confirm the conclusion 
of Nelemans et al. 1 2000) that formation via conservative mass transfer and a common envelope with spiral-in based on energy 
balance or via two such spiral-ins cannot explain the formation of all observed systems. We investigate three different prescrip- 
tions of envelope ejection due to dynamical mass loss with angular-momentum balance and show that they can explain the 
observed masses and orbital periods well. Next, we demand that the age difference of our model is comparable to the observed 
cooling-age difference and show that this puts a strong constraint on the model solutions. However, the scenario in which the 
primary loses its envelope in an isotropic wind and the secondary transfers its envelope, which is then re-emitted isotropically, 
can explain the observed age differences as well. One of these solutions explains the DB-nature of the oldest white dwarf in 
PG 1115+116 along the evolutionary scenario proposed bv lMaxted et alJ l2002all . in which the helium core of the primary 
becomes exposed due to envelope ejection, evolves into a giant phase and loses its hydrogen-rich outer layers. 
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1. Introduction 

Ten double-lined spectroscopic binaries with two white-dwarf 
components are currently known. These binaries have been sys- 
tematically searched for to find possible progenitor systems 
for Type la supernovae, for instance by the SPY (ESO SN la 
Progenitor surveY) project {e.g. Napiwotzki et al. 2001, 2002). 
Short-period double white dwarfs can lose orbital angular mo- 
mentum by emitting gravitational radiation and if the total mass 
of the binary exceeds the Chandrasekhar limit , their eventual 
merge r might produce a supernova of type la Jlben & Tutukovl 

The observed binary systems all have short orbital periods 
that, with one exception, range from an hour and a half to a 
day or two (see Table^, corresponding to orbital separations 
between 0.6 R e and 7 R . The white-dwarf masses of 0.3 M 
or more indicate that their progenitors were (sub)giants with 
radii of a few tens to a few hundred solar radii. This makes 
a significant orbital shrinkage (spiral-in) during the last mass- 
transfer phase necessary and fixes the mechanism for the last 
mass transfer to common-envelope evolution. In such an event 
the envelope of the secondary engulfs the oldest white dwarf 
due to dynamically-unstable mass transfer. Friction then causes 
the two white dwarfs to spiral in towards each other while the 
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envelope is expelled. The orbital energy that is freed due to the 
spiral-in provid es for the necessary energy for the expulsion 
JWebbinkll984l) . 

The first mass transfer phase is usually thought to be either 
another spiral-in or stable and conservative mass transfer. The 
first scenario predicts that the orbit shrinks appreciably during 
the mass transfer whereas the second suggests a widening or- 
bit. Combined with a core-mass -radius relation (e.g. Refsdal 
& Weigert 1970) these scenarios suggest that the mass ratio 
q 2 = M2/M1 of the double white dwarfs is much smaller than 
unity in the first scenario and larger than unity in the second 
scenario. The observed systems all have mass ratio s betw een 
0.70 and 1.28 (Table[T), which led iNelemans et alJ d2000l) to 
conclude that a third prescription is necessary to explain the 
evolution of these systems. They suggested envelope ejection 
due to dynamical mass loss based on angular-momentum bal- 
ance, in which little orbital shrinkage takes place. They used 
analytical approximations to reconstruct the evolution of three 
double white dwarfs and concluded that these three systems 
can only be modelled if this angular-momentum prescription is 
included. 

In this article w e will use the same method as 
INelemans et alJ J2000I) . to see if a stable-mass-transfer episode 
followed by a common envelope with spiral-in can explain the 
observed double white dwarfs. We will improve on their cal- 
culations in several respects. First, we extend the set of ob- 
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served binaries from 3 to 10 systems. Second, we take into ac- 
count progenitor masses for the white dwarf that was formed 
last up to 10 M G and allow them to evo lve beyond core helium 
burning to the asymptotic giant branch. Nelema ns et alJ J2000) 
restricted themselves to progenitor masses of 2.3 M or less 
and did not allow these stars to evolve past the helium flash. 
This was justified because the maximum white-dwarf mass that 
should be created by these progenitors was 0.47 M G , the max- 
imum helium-core mass of a low-mass star and less than the 
minimum mass for a CO white dwarf formed in a spiral-in 
(see Fig.^. The most massive white dwarf in our sample is 
0.71 M Q and cannot have been created by a low-mass star on 
the red-giant branch. Third, we use more sophisticated stellar 
models to reconstruct the evolution of the observed systems. 
This means that the radius of our model stars does not depend 
on the helium-core mass only, but also on total mass of the 
star (see Fig.0. Furthermore, we can calculate the binding en- 
ergy of the hydrogen envelope of our models so that we do not 
need the envelope-structure parameter A enw and can calculate 
the common-envelope parameter a ce directly. Last, because we 
use a full binary-evolution code, we can accurately model the 
stable mass transfer rather than estimate the upper limit for the 
orbital period after such a mass-transfer phase. This places a 
strong constraint on the possible stable-mass-transfer solutions. 
The evolution code also takes into account the fact that the core 
mass of a donor star can grow appreciably during stable mass 
transfer, a fact that alters the relation between the white-dwarf 
mass and the radius of the progenitor mentioned earlier for the 
case of stable mass transfer. 

Our research follows the lines of Nelem ans et alJ |2000), 
calculating the evolution of the systems in reverse order, from 
double white dwarf, via some intermediate system with one 
white dwarf, to the initial ZAMS binary. In Sect.[5]we list the 
observed systems that we try to model. The stellar evolution 
code that we use to calculate stellar models is described in 
Sect.[3] In Sect.|4]we present several grids of single-star mod- 
els from which we will use the helium-core mass, stellar radius 
and envelope binding energy to calculate the evolution during 
a spiral-in. We show a grid of 'basic' models with standard pa- 
rameters and describe the effect of chemical enrichment due to 
accretion and the wind mass loss. We find that these two effects 
may be neglected for our purpose. In Sect.[5]we use the single- 
star models to calculate spiral-in evolution for each observed 
binary and each model star in our grid and thus produce a set 
of progenitor binaries. Many of these systems can be rejected 
based on the values for the common-envelope parameter or or- 
bital period. The remainder is a series of binaries consisting 
of a white dwarf and a giant star that would cause a common 
envelope with spiral-in and produce one of the observed dou- 
ble white dwarfs. In Sect.[6]we model the first mass-transfer 
scenario that produces the systems found in Sect.[5]to com- 
plete the evolution. We consider three possible prescriptions: 
stable and conservative mass transfer, a common envelope with 
spiral-in based on energy balance and envelope ejection based 
on angular-momentum balance. We introduce two variations 
in the latter prescription and show that they can explain the ob- 
served binaries. In addition, we show that the envelope-ejection 
scenario based on angular-momentum balance can also explain 



the second mass-transfer episode. In Sect. 16.41 we include the 
observed age difference in the list of parameters our models 
should explain and find that this places a strong constraint 
on our selection criteria. In Sect.[7]we compare this study to 
earlier work and discuss an alternative formation scenario for 
PG 1 1 15+116. Our conclusions are summed up in Sect.|S] 



2. Observed double white dwarfs 

At present, ten double-lined spectroscopic binaries consisting 
of two white dwarfs have been observed. The orbital periods 
of these systems are well determined. The fact that both com- 
ponents are detected makes it possible to constrain the mass 
ratio of the system from the radial-velocity amplitudes. The 
masses of the components are usually determined by fitting 
white-dwarf atmosphere models to the observed effective tem- 
perature and surface gravity, using mass-radius relations for 
white dwarfs. The values thus obtained are clearly better for the 
brightest white dwarf but less well-constrained than the values 
for the period or mass ratio. It is also harder to estimate the 
errors on the derived mass. In the publications of these obser- 
vations, the brightest white dwarf is usually denoted as 'star 1' 
or 'star A'. Age determinations suggest in most cases that the 
brightest component of these systems is the youngest white 
dwarf. These systems must have evolved through two mass- 
transfer episodes and the brightest white dwarf is likely to have 
formed from the originally less massive component of the ini- 
tial binary (consisting of two ZAMS stars). We will call this 
star the secondary or 'star 2' throughout this paper, whereas the 
primary or 'star 1' is the component that was the initially more 
massive star in the binary. The two components will carry these 
labels throughout their evolution, and therefore white dwarf 1 
will be the oldest and usually the faintest and coldest of the two 
observed components. The properties of the ten double-lined 
white-dwarf systems are listed in Tabled For our calculations 
we will use the parameters that are best determined from the 
Table: P or b, q 2 and M 2 - For M\ we will not use the value listed 
in Tabled but the value Milqi instead. We hereby ignore the 
observational uncertainties in q 2 , because they are small with 
respect to the uncertainties in the mass. In Sects.|5]and|6]we 
will use a typical value of 0.05 M Q JMaxted et alJl2002b¥ for 
the uncertainties in the estimate of the secondary mass. 

Although the cooling-age determinations are strongly de- 
pendent on the cooling model used, the thickness of the hy- 
drogen layer on the surface and the occurrence of shell flashes, 
the cooling-age difference is thought to suffer less from sys- 
tematic errors. The values for At in Ta blefflhave an estimated 
uncertainty of 50% (M axted et alJl2002bl) . The age determina- 
tions of the components of WD 1704+48 la suggest that star 2 
may be the oldest white dwarf, although the age difference 
is small in bo th absolute (20Myr) and relative (^3%) sense 
(Ma xted et alJl2002bi) . Because of this uncertainty we will in- 
troduce an eleventh system with a reversed mass ratio. This 
new system will be referred to as WD 1704+48 lb or 1704b and 
since we assume that the value for M 2 is better determined, we 
will use the following values for this system: M\ = 0.39 M , 
q 2 = 1.43 ± 0.06 and M 2 = q 2 M x = 0.56 M . 
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Table 1. Observed double white dwarfs discussed in this paper. The table shows for each system the orbital period P 01 -b, the 
orbital separation a or b, the masses M\ and M2, the mass ratio q2 = M2/M1, the estimated cooling age of the youngest white 
dwarf T2 and the difference between the cooling ages of the components At. M\ is the mass of the oldest white dwarf and thus 
presumably the original primary. The errors on the peri ods are smaller than the last digit. The values f or a nrh are calculated by th e 
aut hors and mean t to give an indica t ion. R efere nces: (T) | Sa ffer et al Jdl988l). (2)lBergeron et al.Nl989t). ( 3) Braga glia et aljf 1990). 



(4) Marsh ( 1995), (5) Moran et al. ( 1997), (6) Moran et al. ( 1999), (7) Maxted et al. (2000), (8) Bergeron & Liebert (2002), (9) 
iMaxted et al.N2002al) . flOl lMaxted et alJll2002bl) . (I ll lNaniwotzki et alJJ2002h . (12l lKarl et al.N2003al) . (T3l lKarl et alJ (l2003b). 



Note: (a) WD 1704+481a is the close pair of a hierarchical triple. It seems unclear which of the two stars in this pair is the 
youngest (see the text). 



3. The stellar evolution code 

We calculate our models using the STAR S bina r y stellar evo- 
lution code, originally developed by EggletonJ 1 1197 ll 11 972) 
and w ith updated input physics as described in IPols et al. 
Jl995h . Opacity tables are taken from OPAL dlglesias et al. 
1992), compl emented with low-temperature opacities from 
Alexa nder & Ferguson! lfl994). 

The equations for stellar structure and composition are 
solved implicitly and simultaneously, along with an adaptive 
mesh-spacing equation. Because of this, the code is quite sta- 
ble numerically and relatively large timesteps can be taken. As 
a result of the large timesteps and because hydrostatic equilib- 
rium is assumed, the code does not easily pick up short-time- 
scale instabilities such as thermal pulses. We can thus quickly 
evolve our models up the asymptotic giant branch (AGB), with- 
out having to calculate a number of pulses in detail. We thus 
assume that such a model is a good representation of an AGB 
star. 

Convective mixing is modelled by a diffusion equation for 
each of the composition variables, and we assume a mixing- 
length to scale-height ratio l/H p - 2.0. Convective over shoot- 
ing is taken into account as in ISchroder et alJ dl997l) . with 
a parameter <5 ov = 0.12 which corresponds to overshooting 
lengths of about 0.3 pressure scale heights (H p ) and is cali- 
brat ed against accurate stellar data from n on-interacting bina- 
ries JSchroder et allll997HPols et alJll997l) . The code circum- 
vents the helium flash in the degenerate core of a low-mass star 
by replacing the model at which the flash occurs by a model 
with the same total mass and core mass but a non-degenerate 
helium core in which helium was just ignited. The masses of 
the helium and carbon-oxygen cores are defined as the mass 
coordinates where the abundances of hydrogen and helium re- 
spectively become less than 10%. The binding energy of the 
hydrogen envelope of a model is calculated by integrating the 



sum of the internal and gravitational energy over the mass co- 
ordinate, from the helium-core mass M c to the surface of the 
star M s : 




The term t/; nt is the internal energy per unit of mass, that con- 
tains terms such as the thermal energy and recombination en- 
ergy of hydrogen and helium. It has been argued that the bind- 
ing energy of the envelope depends strongly on the definiti on of 
its in ner boundary, i.e. on the definition for M c dDewi & Taurisl 
2000j). This is true for stars with relatively high masses, whose 
cores and envelope binding energies are ill-defined. In our cal- 
culations, however, we deal with relatively low-mass stars on 
the giant branch: these have steep density and composition gra- 
dients at the edge of the core, and as a result the mass of the core 
and binding energy of the envelope in fact are only weakly de- 
pendent on the exact definition of the core. 

We use a version of the code (see Eggleton & Kiseleva- 
Eggleton 2002) that allows for non-conservative binary evolu- 
tion. We use the code to calculate the evolution of both single 
stars and binaries in which both components are calculated in 
full detail. With the adaptive mesh, mass loss by stellar winds 
or by Roche-lobe overflow (RLOF) in a binary is simply ac- 
counted for in the boundary condition for the mass. The spin of 
the stars is neglected in the calculations and the spin-orbit in- 
teraction by tides is switched off. The initial composition of our 
model stars is similar to solar composition: X = 0.70, Y = 0.28 
and Z = 0.02. 

4. Giant branch models 

As we have seen in Sect.^ each of the double white dwarfs that 
are observed today must have formed in a common-envelope 
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event that caused a spiral-in of the two degenerate stars and ex- 
pelled the envelope of the secondary. The intermediate binary 
system that existed before this event, but after the first mass- 
transfer episode, consisted of the first white dwarf (formed 
from the original primary) and a giant-branch star (the sec- 
ondary). This giant is thus the star that caused the common 
envelope and in order to determine the properties of the spiral- 
in that formed each of the observed systems, we need a series of 
giant-branch models. In this section we present a grid of mod- 
els for single stars that evolve from the ZAMS to high up the 
asymptotic giant branch (AGB). For each time step we saved 
the total mass of the star, the radius, the helium-core mass and 
the binding energy of the hydrogen envelope of the star. 

In an attempt to cover all possibilities, we need to take 
into account the effects that can change the quantities men- 
tioned above. We consider the chemical enrichment of the sec- 
ondary by accretion in a first mass-transfer phase and the effect 
a stellar-wind mass loss may have. For each of these changes, 
we compare the results to a grid of 'basic' models with default 
parameters. We keep the overshooting parameter 6 OY constant 
for all these grids, because this effect is unimportant for low- 
mass stars (M < 2.0 M Q ) and its value is well calibrated for 
intermediate-mass stars (see Sect.|3}. 

4.1. Basic models 

In order to find the influence of the effects mentioned above, we 
want to compare the models including these effects to a stan- 
dard. We therefore calculated a grid of stellar models, from the 
zero-age main sequence to high up the asymptotic giant branch 
(AGB), with default values for all parameters. These models 
have solar composition and no wind mass loss. We calculated 
a grid of 199 single-star models with these parameters with 
masses between 0.80 and 10.0 M , with the logarithm of their 
masses evenly distributed. Model stars with masses lower than 
about 2.05 M experience a degenerate core helium flash and 
are at that point replaced by a post-helium-flash model as de- 
scribed in Sect. Because of the large timesteps the code can 
take, the models evolve beyond the point on the AGB where 
the carbon-oxygen core (CO-core) mass has caught up with the 
helium-core mass and the first thermal pulse should occur. 

Figure^ shows the radii of a selection of our grid mod- 
els as a function of their helium-core masses. We used dif- 
ferent line styles to mark different phases in the evolution of 
these stars, depending on their ability to fill their Roche lobes 
or cause a spiral-in and the type of star a common envelope 
would result in. The solid lines show the evolution up the first 
giant branch (FGB), where especially the low-mass stars ex- 
pand much and could cause a common envelope with spiral-in, 
in which a helium white dwarf would be formed. Fig.^i shows 
that low-mass stars briefly contract for core masses around 
0.3 M Q . This is due to the first dredge-up, where the convec- 
tive envelope deepens down to just above the hydrogen burning 
shell and increases the hydrogen abundance there. The contrac- 
tion happens when the hydrogen-burning shell catches up with 
this composition discontinuity. After ignition of helium in the 
core, all stars shrink and during core helium burning and the 
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Fig. 1. Core-mass -radius relations for the 'basic' models, as 
described in the text. The lines show the logarithm of the ra- 
dius of the stars as a function of the helium-core mass. Upper 
panel (a): low-mass grid models with masses of 0.91, 1.01, 
1.14, 1.30, 1.48, 1.63, 1.81 and 2.00 M Q . Lower panel (b): high- 
mass grid models with masses of 2.00, 2.46, 2.79, 3.17, 3.60, 
4.09, 4.65, 5.28 and 6.00 M@. The 2 M model is plotted in both 
panels throughout as a solid line for easier comparison. The 
other models are shown as solid lines on the first giant branch 
(FGB), where they could cause a common envelope resulting in 
a spiral-in and creating a helium white dwarf. The dashed lines 
show the asymptotic giant branch (AGB), where a common en- 
velope would lead to the formation of a carbon-oxygen white 
dwarf. Dotted lines are parts of the evolution where the stars 
either are smaller than at the tip of the FGB (at lower radii) or 
where their envelope binding energies become positive on the 
AGB (at large radii). 



first phase of helium fusion in a shell, their radii are smaller 
than at the tip of the FGB. This means that these stars could 
never start filling their Roche lobes in this stage. These parts of 
the evolution are plotted with dotted lines. Once a CO core is 
established, the stars evolve up the AGB and eventually get a 
radius that is larger than that on the FGB. The stars are now ca- 
pable of filling their Roche lobes again and cause a common en- 
velope with spiral-in. In such a case we assume that the whole 
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Fig. 2. Helium-core-mass -orbital period tracks for the 'basic' 
models. The lines show the logarithm of the orbital period at 
which the Roche lobe is filled for grid models with masses of 
1.01, 1.27, 1.59, 2.00, 2.52, 3.17, 3.99, 5.02 and 6.32M Q . The 
period was obtained from the radius of the model star, under 
the assumption that it fills its Roche lobe and the companion 
has a mass equal to the helium-core mass of the model. This 
way, the system would undergo a spiral-in that would lead to 
a binary with mass ratio q — 1 . The line styles have the same 
meaning as in Fig.^ 



helium core survives the spiral-in and that the helium burning 
shell will convert most of the helium to carbon and oxygen, 
eventually resulting in a CO white dwarf, probably with an at- 
mosphere that consists of a mixture of hydrogen and helium. 
This part of the evolution is marked with dashed lines. Fig.^ 
shows that the most massive models in our grid have a decreas- 
ing helium-core mass at some point on the AGB. This happens 
at the so-called second dredge-up, where the convective mantle 
extends inward, into the helium core and mixes some of the he- 
lium from the core into the mantle, thereby reducing the mass 
of the core. Models with masses between about 1.2 and 5.6 M 
expand to such large radii that the binding energy of their hy- 
drogen envelopes become positive. In Sect.|5]we are looking for 
models that can cause a spiral-in based on energy balance in the 
second mass-transfer phase, for which purpose we require stars 
that have hydrogen envelopes with a negative binding energy. 
A positive binding energy means that there is no orbital energy 
needed for the expulsion of the envelope and thus the orbit will 
not shrink during a common envelope caused by such a star. We 
have hereby implicitly assumed that the recombination energy 
is available during common-envelope ejection. 

To give some idea what kind of binaries can cause a spiral- 
in and could be the progenitors of the observed double white 
dwarfs, we converted the radii of the stars displayed in Fig.^ 
into orbital periods of the pre-common-envelope systems. To 
do this, we assumed that the Roche-lobe radius is equal to the 
radius of the model star, and that the mass of the companion is 
equal to the mass of the helium core of the model. This is jus- 
tified by Tabled where the geometric mean of the mass ratios 
is equal to 1.03. The result is shown in Fig.|2] 
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Fig. 3. The logarithm of the binding energy of the 'basic' 
model stars as a function of the helium-core mass. The grid 
models with masses of 0.91, 1.01, 1.14, 1.30, 1.48, 1.63, 1.81, 
2.00, 2.46, 2.79, 3.17, 3.70, 4.09, 4.65, 5.28, 6.00 and 6.82 M Q 
are shown. The 2.00 M Q model is drawn as a solid line, the line 
styles for the other models have the same meaning as in Fig.^ 
The parts where the envelope binding energy is zero (before a 
helium core develops) or positive are not shown. 

In Sect.|5]we will need the efficiency parameter a ce of each 
common-envelope model to judge whether that model is ac- 
ceptable or not. In order to calculate this parameter we must 
know the binding energy of the hydrogen envelope of the pro- 
genitor star (see Eq.0}, that is provided by the evolution code 
as shown in Eq.^ The envelope binding energy is therefore an 
important parameter and we show it for a selection of models 
in Fig.|3] again as a function of the helium-core mass. Because 
the binding energy is usually negative, we plot the logarithm of 
-f7b, e - The phases where the envelope binding energy is non- 
negative are irrelevant for our calculations of ff ce and therefore 
not shown in the Figure. 

Many common-envelope calculations in the literature use 
the so-called envelope-structure parameter /l env to estimate the 
envelope binding energy from basic stellar parameters in case 
a detailed model is not available 



De Kool et al. (1987) suggest that A ew/ » 0.5. Since we cal- 
culate the binding energy of the stellar envelope, we can in- 
vert Eq.|3and calculate A em (because we mostly consider low- 
mass stars, the binding energy and hence A em does not depend 
strongly on the definition of the core mass, see the discussion 
below Eq.^and see also Dewi & Tauris 2000). Figure@]shows 
the results of these calculations as a function of the helium core 
mass, for the same selection of models as in Fig. [5] We see that 
a value of /l env = 0.5 is a good approximation for the lower 
FGB of a low-mass star, or the FGB of a higher-mass star. A 
low-mass star near the tip of the first giant branch has a struc- 
ture parameter between 0.5 and 1.5 and for most stars A em in- 
creases to more than unity rather quickly, especially when the 
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Fig. 4. The envelope-structure parameter A em for the 'basic' 
models, as a function of the helium-core mass. The same grid 
models are shown as in Fig. [2] The meaning of the line styles 
is explained in the caption of Fig.^ 



stars expand to large radii and the binding energies come close 
to zero. 

4.2. Chemical enrichment by accretion 

The secondary that causes the common envelope may have 
gained mass by accretion during the first mass-transfer phase. 
If this mass transfer was stable, the secondary has probably 
accreted much of the envelope of the primary star. The deep- 
est layers of the envelope of the donor are usually enriched 
with nuclear burning products, brought up from the core by a 
dredge-up process. This way, the secondary may have been en- 
riched with especially helium which, in sufficiently large quan- 
tities, can have an appreciable effect on the opacity in the en- 
velope of the star and thus its radius. This would change the 
core-mass -radius relation of the star and the common enve- 
lope it causes. 

To see whether this effect is significant, we considered a 
number of binary models that evolved through stable mass 
transfer to produce a white dwarf and a main-sequence sec- 
ondary. The latter had a mass between 2 and 5 M in the cases 
considered, of which 50-60% was accreted. We then took this 
secondary out of the binary and let it evolve up the asymptotic 
giant branch, to the point where the code picks up a shell in- 
stability and terminates. We then compared this final model to 
a model of a single star with the same mass, but with solar 
composition, that was evolved to the same stage. In all cases 
the core-mass -radius relations coincide with those in Fig.^ 
When we compared the surface helium abundances of these 
models, after one or two dredge-ups, we found that although 
the abundances were enhanced appreciably since the ZAMS, 
they were enhanced with approximately the same amount and 
the relative difference of the helium abundance at the surface 
between the different models was always less than 1.5%. In 
some cases the model that had accreted from a companion had 
the lower surface helium abundance. 



tip of FGB 

+ tip of AGB 




Mj (M ) 

Fig. 5. The fraction of mass lost at two moments in the evolu- 
tion of a star as a function of its initial mass, for the three dif- 
ferent wind strengths (C rm i = 0.2, 0.5 and 1.0) used in the grid. 
This fraction is shown for the tip of the FGB (dashed lines and 
crosses), and the 'tip of the AGB' (dotted lines and plusses). 
See the text for details. 

The small amount of helium enrichment due to accretion 
gives rise to such small changes in the core mass -radius re- 
lation, that we conclude that this effect can be ignored in our 
common-envelope calculations in Sect.[5] 

4.3. Wind mass loss 

The mass loss of a star by stellar wind can change the mass of 
a star appreciably before the onset of Roche-lobe overflow, and 
the mass loss can influence the relation between the core mass 
and the radius of a star. From Fig.[T]it is already clear that this 
relation depends on the total mass of the star. In this section, we 
would therefore like to find out whether a conservative model 
star of a certain total mass and core mass has the same radius 
and envelope binding energy as a model with the same total 
mass and core mass, but that started out as a more massive star, 
has a strong stellar wind and just passes by this mass on its 
evolution down to even lower masses. We calculated a small 
grid of models with ten different initial masses between 1.0 M 
and 8.0 M , ev enly spread in logM and included a Reimers 
type mass loss jReimersI 19751) of variable strength: 



(3) 



where we have used the values C rm i = 0.2, 0.5 and 1.0. The 
basic models of Sect. 14.11 are conservative and therefore have 
Crmi = 0. The effect of these winds on the total mass of the 
model stars in our grid is displayed in Fig.|5] It shows the frac- 
tion of mass lost at the tip of the first giant branch (FGB) and 
the 'tip of the asymptotic giant branch' (AGB). The first mo- 
ment is defined as the point where the star reaches its largest 
radius before helium ignites in the core, the second as the point 
where the radius of the star reaches its maximum value while 
the envelope binding energy is still negative. Values for both 
moments are plotted in Fig.|5]for each non-zero value of C rm i in 
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Fig. 6. Comparison of a selection from the small grid of models 
with a stellar wind. The models displayed have masses of 1 .0, 
1.6, 2.5, 3.2, 4.0, 5.0 and 6.3 M Q . The wind strength parameters 
are C rm i = 0.0 (dotted lines) and C rm i = 1.0 (solid lines, the 
strongest mass loss in the grid). Stars with mass loss are usually 
larger, but for models of 4.0 M or more this effect becomes 
negligible. The 1.0 M G model loses so much mass that it never 
ignites helium in the core. 

the grid. For the two models with the lowest masses the highest 
mass-loss rates are so high that the total mass is reduced suf- 
ficiently on the FGB to keep the star from igniting helium in 
the core, and the lines in the plot coincide. Stars more massive 
than 2 M G have negligible mass loss on the FGB, because they 
have non-degenerate helium cores so that they do not ascend 
the FGB as far as stars of lower mass. Their radii and luminosi- 
ties stay relatively small, so that Eq.[3] gives a low mass loss 
rate. For stars of 4 M or more, the mass loss is diminutive and 
happens only shortly before the envelope binding energy be- 
comes positive. We can conclude that for these stars the wind 
mass loss has little effect on the core-mass -radius relation. 

The core-mass -radius relations for a selection of the mod- 
els from our wind grid are shown in Fig. [6] The Figure com- 
pares models without stellar wind with models that have the 
strongest stellar wind in our grid (C rm i = 1.0). Models with the 
other wind strengths would lie between those shown, but are 
not plotted for clarity. The greatest difference in Fig.[6]is in the 
1.0 M model. The heavy mass loss reduces the total mass of 
the star to 0.49 M e on the first giant branch, so that the star is 
not massive enough to ignite helium in the core. Fig.[6] shows 
that models with mass loss are larger than conservative mod- 
els for the same core mass, as one would expect from Fig.^ 
This becomes clear on the FGB for stars that have degenerate 
helium cores, because they have large radii and luminosities 
and lose large amounts of mass there. For stars more massive 
than about 2 M the mass loss becomes noticeable on the AGB. 
Stars of 4 M or more show little difference in Fig. [6] The enve- 
lope binding energies have similar differences in the same mass 
regions. 

The question is whether the properties of the model with 
reduced mass due to the wind are the same as those for a con- 
servative model of that mass. In order to answer this question, 



we have compared the models from the 'wind grid' to the ba- 
sic, conservative models. As the wind reduces the total mass 
of a model star, it usually reaches masses that are equal to that 
of several models in the conservative grid. As this happens, we 
interpolate linearly within the mass-losing model to find the 
exact moment where its mass equals the mass of the conser- 
vative model. We then use the helium-core mass of the inter- 
polated mass-losing model to find the moment where the con- 
servative model has the same core mass and we calculate its 
radius and envelope binding energy, again by linear interpola- 
tion. This way we can compare the two models at the moment 
in evolution where they have the same total mass and the same 
core mass. This comparison is done in Fig.0 Figure^ directly 
compares the radii of the two sets of models, in Fig.03 the ratio 
of the two radii is shown. 

Of the data points in Fig. [TJ) 83% lie between 0.9 and 1.1 
and 61% between 0.95 and 1.05. For the wind models with 
Crmi = 0.2 these numbers are 99% and 97%, and for the mod- 
els with Crmi = 0.5 they are 94% and 85% respectively. As can 
be expected, the models that have a lower - and perhaps a more 
realistic - mass-loss rate compare better to the conservative 
models. We see in Fig.Et that many of the points that lie farther 
from unity need only a small shift in core mass to give a perfect 
match. This shift is certainly less than 0.05 M , which is what 
we will adopt for the uncertainty of the white-dwarf masses 
in Sect.|5] We conclude here that there is sufficient agreement 
between a model that reaches a certain total mass because it 
suffers from mass loss and a conservative model of the same 
mass. The agreement is particularly good for stars high up on 
the FGB or AGB, where the density contrast between core and 
envelope is very large. 

5. Second mass-transfer phase 

For the formation of two white dwarfs in a close binary sys- 
tem, two phases of mass transfer must happen. We will call 
the binary system before the first mass transfer the initial bi- 
nary, with masses and orbital period Mn, and P\, If one 
considers mass loss due to stellar wind before the first mass- 
transfer episode, these parameters are not necessarily equal to 
the ZAMS parameters, especially for large 'initial' periods. 
The binary between the two mass-transfer phases is referred 
to as the intermediate binary with M\ m , Mi m and P m . After the 
two mass-transfer episodes, we obtain the final binary with pa- 
rameters M[f, M%[ and Pf, that should correspond to the values 
that are now observed and listed in Tabled The subscripts ' 1 ' 
and '2' are used for the initial primary and secondary as defined 
in Sect.El 

In the first mass transfer, the primary star fills its Roche lobe 
and loses mass, that may or may not be accreted by the sec- 
ondary. This leads to the formation of the intermediate binary, 
that consists of the first white dwarf and a secondary of un- 
known mass. In the second mass-transfer phase, the secondary 
fills its Roche lobe and loses its envelope. The second mass 
transfer results in the observed double white dwarf binaries that 
are listed in Tabled and must account for significant orbital 
shrinkage. This is because the youngest white dwarf must have 
been the core of its progenitor, the secondary in the interme- 
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Fig. 7. Comparison of a selection of grid models with C lm i = 
1.0 with initial masses of 1.3, 1.6, 2.0, 2.5 and 3.2 M to the 
basic models (C rm i = 0.0). Upper panel (a): Comparison of 
the radius of the models with a stellar wind (solid lines) and 
the radius of a basic model with the same mass and core mass 
(plusses). Lower panel (b): The fraction of the radius of the 
wind model R w over the radius of the basic model R\, with 
the same total and core mass. Each data point corresponds 
to a point in the upper panel. Of the data points in the upper 
panel, 7 out of 143 (5%) lie outside the plot boundaries in the 
lower panel. The dashed horizontal lines show the region where 
agreement is better than 10%, where 83% of the data points lie. 
The 1 .0 M model was left out because there are only a few ba- 
sic models with lower mass, the higher-mass models were left 
out because they lose very little mass (see Fig.[5}- 



diate binary. Stars with cores between 0.3 and 0.7 M usually 
have radii of several tens to several hundreds of solar radii, and 
the orbital separation of the binaries they reside in must be even 
larger than that. The orbital separation of the observed systems 
is typically only in the order of a few solar radii (Table[Q. Giant 
stars with large radii have deep convective envelopes and when 
such a star fills its Roche lobe, the ensuing mass transfer will 
be unstable and occur on a very short, dynamical timescale, es- 
pecially if the donor is much more massive than its companion. 
It is thought that the envelope of such a star can engulf its com- 



panion and this event is referred to as a common envelope. The 
companion and the core of the donor orbit inside the common 
envelope and drag forces will release energy from the orbit, 
causing the orbit to shrink and the two degenerate stars to spiral 
in. The freed orbital energy will heat the envelope and eventu- 
ally expel it. This way, the hypothesis of the common envelope 
with spiral-in can phenomenologically explain the formation of 
close double-white-dwarf binaries. 

5.1. The treatment of a spiral-in 

In order to estimate the orbital separation of the post-common 
envelope system quantitatively, it is often assumed that the or- 
bital energy of the system is decreased by an amount that is 
equal to t he binding energy of the envelope of the donor star 
(Webb inkll98l : 



t/b,e 



GM lf M 2i GM lm M 2n 



2flf 



2a„ 



(4) 



The parameter a ce is the common-envelope parameter that ex- 
presses the efficiency by which the orbital energy is deposited 
in the envelope. Intuitively one would expect that a ce « 1. 
However, part of the liberated orbital energy might be radiated 
away from the envelope during the process, without contribut- 
ing to its expulsion, thereby lowering a ce . Conversely, if the 
common-envelope phase would last long enough that the donor 
star can produce a significant amount of energy by nuclear fu- 
sion, or if energy is released by accretion on to the white dwarf, 
this energy will support the expulsion and thus increase a ce . 

In the forward calculation of a spiral-in the final orbital 
separation Of depends strongly on the parameter ar ce , which 
must therefore be known. In this section we will try to estab- 
lish the binary systems that were the possible progenitors of 
the observed double white dwarfs and we will therefore per- 
form backward calculations. The advantage of this is that we 
start as close as possible to the observations thus introducing 
as little uncertainty as possible. The problem with this strategy 
is that we do not know the mass of the secondary progenitor 
beforehand. We will have to consider this mass as a free pa- 
rameter and assume a range of possible values for it. The grid 
of single-star models of Sect.0]provides us with the total mass, 
core mass, radius and envelope binding energy at every mo- 
ment of evolution, for a range of total masses between 0.8 and 
10 M e . It is then not necessary to know the common-envelope 
parameter, and we can in fact use the binding energy to cal- 
culate the value for ff C e that is needed to shrink the orbit of a 
model to the period of the observed double white dwarf. The 
accuracy with which a ce can be calculated scales directly with 
the accuracy of the binding energy. As discussed below Eq.^ 
[/b.e and as a consequence a ce: do not depend strongly on the 
exact way in which the core is defined, for the low-mass stars 
which determine the outcome of our calculations (see Table|5] 
below). We make two assumptions about the evolution of the 
two stars during the common envelope to perform these back- 
ward calculations: 

1 . the core mass of the donor does not change, 

2. the mass of the companion does not change. 
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The first assumption will be valid if the timescale on which the 
common envelope takes place is much shorter than the nuclear- 
evolution timescale of the giant donor. This is certainly true, 
since the mass transfer occurs on the dynamical timescale of 
the donor. The second assumption is supported firstly by the 
fact that the companion is a white dwarf, a degenerate ob- 
ject that has a low Eddington accretion limit and is further- 
more difficult to hit directly by a mass stream from the donor. 
The whit e dwarf could accr ete matter in the Bondi-Hoyle fash- 
ion JBondi & Hovlelfl944h . This would not change the mass 
of the white dwarf significantly but could release appreciable 
amounts of energy. Secondly, a common envelope is estab- 
lished very shortly after the beginning of the mass transfer, so 
that the mass stream disappears and the white dwarf is orbiting 
inside the fast-expanding envelope rather than accreting mass 
from the donor. In the terminology used here, the second as- 
sumption can be written asMi m = M\{ . 

From the two assumptions above it follows that the mass of 
the second white dwarf, the one that is formed in the spiral-in, 
is equal to the helium-core mass of the donor at the moment 
it fills its Roche lobe. There is therefore a unique moment in 
the evolution of a given model star at which it could cause a 
common envelope with spiral-in and produce a white dwarf of 
the proper mass. Recall from Fig.^i that although the second 
dredge-up reduces notably the helium-core mass of the more 
massive models in the grid, there is no overlap in core mass in 
the phases where the star could fill its Roche lobe on the first 
giant branch (solid lines) or asymptotic giant branch (dashed 
lines). The moment where the model star could produce a white 
dwarf of the desired mass in a common envelope with spiral-in 
is therefore defined by two conditions: 

1 . the helium-core mass of the model reaches the mass of the 
white dwarf, 

2. the model star has its largest radius so far in its evolution. 

The second restriction is necessary because stars can shrink ap- 
preciably during their evolution, as noted in Sect. 14.11 If the 
core of a model star obtains the desired mass at a point in the 
evolution where the star is smaller than it has been at some 
point in the past, then the star cannot fill its Roche lobe at the 
right moment to produce a white dwarf of the proper mass and 
therefore this star cannot be the progenitor of the white dwarf. 
This way, each model star has at most one moment in its evo- 
lution where it could fill its Roche lobe and produce the ob- 
served double white dwarf. If such a moment does not exist, 
the model star is rejected as a possible progenitor of the second 
white dwarf. 

If the model star could be the progenitor of the youngest 
white dwarf in the observed system, the computer model gives 
us the radius of the donor star, that must be equal to the Roche- 
lobe radius. Under the assumption that the mass of the first 
white dwarf does not change in the common envelope, the mass 
ratio of the two stars q 2m = M 2m /Mi m and the Roche-lobe ra- 
dius of the secondary star RR\ 2m give us the orbital sep aration 
before the spiral-in a m , where we use the fit by Eggleton ( 1983) 



^R12m - An 



0.49 <?£ 



0.6^ + ln(l + ^ 3 ) 



< q 2m < oo. 



(5) 



Kepler's law finally provides us with the orbital period P m of 
the intermediate system. The stellar model also gives the bind- 
ing energy of the envelope of the donor t/b, e at the onset of 
the common envelope and we can use Eq.0]to determine the 
common-envelope parameter a ce . We will use a ce to judge the 
validity of the model star to be the progenitor of the second 
white dwarf. There are several reasons why a numerical so- 
lution can be rejected. Firstly, the proposed donor could be a 
massive star with a relatively small radius. Then a m will be 
small and it might happen that a m < Of j^-, so that a ce < 0. 
This means that energy is needed to change the orbit from a m 
to «f, or even that a m < Of and a spiral-in (if it can be called 
that) to the desired orbit will not lead to expulsion of the com- 
mon envelope. Secondly, as mentioned above, a ce is expected 
to be close, though not necessarily equal, to unity. However if 
the parameter is either much smaller or much larger than 1, 
we will consider the spiral-in to be 'physically unbelievable' . 
We arbitrarily chose the boundaries between which or ce must 
lie for a believable spiral-in to be a factor of ten either way: 
0.1 < a ce < 10. We think that the actual value for a ce should 
be more constrained than that because common-envelope evo- 
lution is thought to last only a short time so that there is little 
time to generate or radiate large amounts of energy, but keep 
the range as broad as it is to be certain that all possible progen- 
itor systems are considered in our sample. 



5.2. Results of the spiral-in calculations 

We will now apply the stellar models of Sect. l4.ll as described 
in the previous section to calculate potential progenitors to the 
observed double white dwarfs as listed in Tabled As input pa- 
rameters we took the values Pf = P OI \, and M 2 f = M 2 from 
the table, and assumed that Mu = M 2 /q 2 , where q 2 is the ob- 
served mass ratio listed in Tabled We thus ignore for the mo- 
ment any uncertainty in the observed masses. Figure[8] shows 
the orbital period P m as a function of the secondary mass M 2m . 
Each symbol is a solution to the spiral-in calculations and rep- 
resents an intermediate binary system that consists of the first 
white dwarf of mass M\ m - M\f, a companion of mass M 2m 
and an orbital period P m . The secondary of this system will 
fill its Roche lobe at the moment when its helium-core mass 
is equal to the mass of the observed white dwarf M 2 f, and can 
thus produce the observed double-white-dwarf system with a 
common-envelope parameter that lies between 0. 1 and 10. 

The solutions for each system in Fig.[8]seem to lie on curves 
that roughly run from long orbital periods for low-mass donors 
to short periods for higher-mass secondaries. This is to be ex- 
pected, partially because higher-mass stars have smaller radii 
at a certain core mass than stars of lower mass (see Fig.[0 and 
thus fill their Roche lobes at shorter orbital periods, but mainly 
because the orbital period of a Roche-lobe filling star falls off 
approximately with the square root of its mass. The Figure also 
shows gaps between the solutions, for instance for WD 0957- 
666 and WD 1704+48 la, between progenitor masses M 2m of 
about 2 and 2.5 M Q . These gaps arise because the low-mass 
donors on the left side of the gap ignite helium degenerately 
when the core mass is 0.47 M , after which the star shrinks, 
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+0135 ♦Q136 O0957 

k 1101 □ 11 15 A1204 

»1349 il414 * 1704a 
* 1704b .2209 



M 2m < M o) 

Fig. 8. Results of the spiral-in calculations, each individual 
symbol is a solution of the calculations and thus represents 
one pre-CE binary. The figure shows the logarithm of the or- 
bital period of the intermediate binary P m as a function of the 
secondary mass Mim- Different symbols represent different ob- 
served systems, as explained in the legend. System 1704a is 
the system listed in Tabled 1704b is the same system, but with 
the reverse mass ratio. For solutions with Mim < 2.5 M , only 
every third solution is plotted for clarity. Around M 2m = 1.2 
and log P m = 2.8 the symbols of WD 0135-052, WD 0136+768 
and WD 1204+450 overlap due to the fact that they have sim- 
ilar white-dwarf masses. For comparison we show the lines of 
the solutions for (top to bottom) WD 01 36+768 WD 09 57-666 
and WD 1 101+364 taken from iNelemans et alJ ll2000l) . as de- 
scribed in the text. 



whereas for stars with masses close to 2 M helium ignition 
is non-degenerate and occurs at lower core masses, reaching a 
minimum for stars with a mass of 2.05 M , where helium ig- 
nition occurs when the helium-core mass amounts to 0.33 M 
(see Fig.Q. Thus, for white dwarfs with masses between 0.33 
and 0.47 M Q there is a range of masses between about 1.5 and 
3 M Q for which the progenitor has just ignited helium in the 
core, and thus shrunk, when it reaches the desired helium-core 
mass. 

The dip and gap in Fig.|8]for WD 1101+364 (with M 2f « 
O.29M ) around M 2m = 1 .8 M can be attributed to the first 
dredge-up that occurs for low-mass stars (M < 2.2 M ) early 
on the first giant branch. Stars with these low masses shrink 
slightly due to this dredge-up that occurs at core masses be- 
tween about 0.2 and 0.33 M , the higher core masses for the 
more massive stars (see Fig.[^). Stars at the low-mass (M2 m ) 
side of the gap obtain the desired core mass just after the 
dredge-up, are relatively small and fill their Roche lobes at 
short periods. Stars with masses that lie in the gap reach that 
core mass while shrinking and cannot fill their Roche lobes 
for that reason. Stars at the high-mass end of the gap fill their 
Roche lobes just before the dredge-up so that this happens 
when they are relatively large and therefore this happens at 
longer orbital periods. 

For comparison we display as solid lines in Fig. [8] 
the results for the systems WD 0136+768, WD 0957-666 



+0135 ♦0136 O0957 

x 1101 =1115 M304 

• 1349 M414 * 1704a 

* 1704b - 2209 



M 2m (Mo) 

Fig. 9. The logarithm of the common-envelope parameter a ce 
for the solutions of the spiral-in calculations shown in Fig. [8] 
Different symbols represent different observed systems. For 
M2m < 2.5 M every third solution is plotted only. 



and WD 1101+364 ( from top to bottom), as found by 
Nelemans et al and shown in their Fig. 1. The differ- 

ences between their and our results stem in part from the fact 
that the values for the observed masses have been updated by 
observations since their paper was published. To compensate 
for this we include dashed lines for the two systems for which 
this is the case. The dashed lines were calculated with their 
method but the values for the observed masses as listed in this 
paper. By comparing the lines to the symbols for the same sys- 
tems, we see that they lie in the same region of the plot and 
in the first order approach they give about the same results. 
However, the slopes in the two sets of results are clearl y dif- 
ferent. This can be attributed to the fact that INelemans et all 
(2000) used a power law to describe the radius of a star as a 
function of its core mass only. The change in orbital period 
with mass in their calculations is the result of changing the total 
mass in Kepler's law. Furthermore, they assumed that all stars 
with masses between 0.8 and 2.3 M have a solution, whereas 
we find limits and gaps, partially due to the fact that we take 
into account the fact that stars shrink and partially because in 
Fig.|8]only solutions with a restricted a ce are allowed. On the 
other hand, we allow stars more massive than 2.3 M Q as possi- 
ble progenitors. 

In Fig.|9] we display the common-envelope parameter a ce 
for a selection of the solutions with 0.1 < a ce < 10 as a func- 
tion of the unknown intermediate secondary mass M2 m . Each 
of the plot symbols has a corresponding symbol in Fig.|S] To 
produce these two figures, we have so far implicitly assumed 
that the masses of the two components are exact, so that there 
is at most one acceptable solution for each progenitor mass. 
This is of course unrealistic and it might keep us from find- 
ing an acceptable solution. At this stage we therefore introduce 
an uncertainty on the values for Mi in Tableland take M 2 f = 
M 2 - 0.05 M , M 2 - 0.04 M Q , . . . , M 2 + 0.05 M . Meanwhile 
we assume that the mass ratio and orbital period have negligi- 
ble observational error, because these errors are much smaller 
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M 2m ( M 0) 

Fig. 10. Results of the spiral-in calculations. This figure is sim- 
ilar to Fig.|8]and shows the orbital period of the pre-CE system 
as a function of the secondary mass. The solutions for each 
system are plotted in a separate panel, as labelled in the upper- 
right corner. All solutions with acceptable a ce: are plotted. The 
number of solutions for each system is shown in the lower-right 
corner. We assumed an uncertainty in M 2 t of 0.05 M and cal- 
culated M lt - using q 2 . 



than those on the masses, and obtain the mass for the first white 
dwarf from M\f = M 2 f/q 2 . Thus we have 1 1 pairs of values for 
M\f and M 2 f for each observed system, which we use as input 
for our spiral-in calculations. The results are shown in Fig.llQI 

If we compare Fig. [3] and Fig.EH we see that the wider 
range in input masses results in a wider range of solutions, 
similar to those we found in Fig. [8] but extended in orbital pe- 
riod. This can be understood qualitatively, since lowering the 
white-dwarf mass demands a lower helium-core mass in the 
progenitor and thus a less evolved progenitor with a smaller ra- 
dius at the onset of Roche-lobe overflow. Conversely, higher 
white-dwarf masses need more evolved progenitors that fill 
their Roche lobes at longer orbital periods. The introduction 
of this uncertainty clearly results in a larger and more realis- 
tic set of solutions for the spiral-in calculations and therefore 
should be taken into account. 

Each system in Fig.^|is a possible progenitor of one of 
the ten observed double white dwarfs listed in Tabled We now 
turn to the question whether and how these intermediate sys- 
tems can be produced. 

6. First mass-transfer phase 

The solutions of the spiral-in calculations we found in the pre- 
vious section are in our nomenclature intermediate binaries, 
that consist of one white dwarf and a non-degenerate compan- 
ion. In this section we will look for an initial binary that con- 
sists of two zero-age main-sequence (ZAMS) stars of which 
the primary evolves, fills its Roche lobe, loses its hydrogen en- 
velope, possibly transfers it to the secondary, so that one of 
the intermediate binaries of Fig.^3 1S produced. The nature 
of this first mass transfer is a priori unknown. In the follow- 



ing subsections we will consider (1) stable and conservative 
mass transfer that will result in expansion of the orbit in most 
cases, (2) a common envelope with spiral-in based on energy 
balance (see Eq.@} that usually gives rise to appreciable orbital 
shrinkage and (3) envelope ejection due to dynamically unsta- 
ble mass loss based on angular-momentu m balance, as intro- 
duced by lPaczvri ski & Z iolkowskil i 19671) and already used by 
Nele mans et alJ (|2000) for the same purpose, which can take 
place without much change in the orbital period. 

6.1. Conservative mass transfer 

In this section we will find out which of the spiral-in solu- 
tions of Fig.^|may be produced by stable, conservative mass 
transfer. We use the binary evolution code described in Sect.|3| 
For simplicity, we ignore stellar wind and the effect of stellar 
spin on the structure of the star. Because we assume conserva- 
tive evolution, the total mass of the binary is constant, so that 
M\{ + Mjx — M\ m + M 2m , where the last two quantities are 
known. Also, we ignore angular momentum exchange between 
spin and orbit by tidal forces, so that the orbital angular mo- 
mentum is conserved. This implies that 

= I M u M 2 i ? 
Pi \M lm M 2 J ' 

Because of the large number of possible intermediate sys- 
tems we will first remove all such systems for which it can a 
priori be shown that they cannot be produced by conservative 
mass transfer. These systems have orbital periods that are either 
too short or too long to be formed this way. We can find a lower 
limit to the intermediate period as a function of secondary mass 
M 2m using the fact that the total mass of the initial system must 
be equal to the sum of the mass of the observed white dwarf 
M\ and M 2m . We distributed this mass equally over two ZAMS 
stars and set the Roche-lobe radii equal to the two ZAMS radii. 
By substituting the initial and desired masses in Eq.[6]we find a 
lower limit to the period of the intermediate binary, which we 
will call P m j n . 

For each mass of the intermediate binary, an upper limit to 
the intermediate period P m can be estimated in two steps, as 
follows. First, we follow Nelem ans et ail (2000) in noting that 
a maximum period after mass transfer is reached for an initial 
mass ratio <72i,opt = 0.66. 1 To reach the highest value for P m 
the initial mass ratio must be optimal as defined above, and 
in addition the initial period must be the longest period P; i0p t 
for which the mass transfer is stable. We use the conditions by 
iHurlev etal] f2000) who define this point as the moment where 
the mass of the convective envelope Mqe exceeds a certain frac- 
tion of the total mass of the hydrogen envelope Me for the first 
time: 



M CE 
M CE 



■Me 



Mnnpt < 1.995 M 0> 



\ Me, M h , 0£L > 1-995 M, 



(7) 



1 Nelemans et al. (2000) find ^2i,opt = 0.62, because the equation 
they use for the Roche lobe is different from the Eq.|5|that we use. 
They also use a different condition for the stability of mass transfer 
from an evolving star. Our calculations show that if we would use 
Eq.|7| instead, the upper limit to P m would drop, so that we can be 
confident that the limit we find is indeed an upper limit. 
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1 2 3 4 5 

Fig. 11. Results of the spiral-in calculations for WD 0957-666 
with period limits for a conservative first mass transfer. This 
figure contains the same data as the third panel in Fig.llOKsvm- 
bols) plus the period limits P m ; n and P max (dashed lines). The 
solutions that lie between these limits are shown in black, the 
others in grey. See the main text for details. 

for Z = 0.02. We then interpolate in our grid of stellar mod- 
els of Sect.|4]to find the radius of a star with the desired mass 
at the base of the giant branch (Rbgb)- By assuming that this 
radius is equal to the Roche-lobe radius and using Eq.|5] we 
find the optimum initial period Pi j0pt . For a given binary mass, 
a unique initial binary system is thus defined by opt , M2i t0pt 
and Piopt- We then use Eq.|6]to find for this initial system the 
maximum intermediate period, which we will call f max . All in- 
termediate systems, that result from our spiral-in calculations 
and have longer orbital periods than f max , can not be formed 
by conservative mass transfer. 

The lower and upper limits for the orbital period between 
which a conservative solution must lie for WD 0957-666 are 
shown in Fig.^Jtogether with the intermediate systems found 
from the spiral-in calculations. Black dots represent solutions 
that lie between the limits and could match the outcome of a 
conservative model, grey dots lie outside these limits and can- 
not be created by conservative mass transfer. There is a slight 
difference between the dashed lines and the division between 
filled and open symbols in the Figure, because the spiral-in 
solutions are shown with the uncertainty in the masses de- 
scribed in the previous section, whereas the period limits are 
only shown for the measured M2 and q-i (see Tabled for clar- 
ity. 

After selecting the spiral-in solutions that lie between these 
period limits for all eleven systems, we find that such solu- 
tions exist for only six of the observed binaries, as shown in 
Fig-El We tried to model these intermediate systems with the 
binary evolution code described in Sect.[3] Because of the large 
number of allowed spiral-in solutions for WD 0957-666 and 
WD 1 101+364, we decided to model about half of the solutions 
for these two systems and all of the solutions for the other four. 
Because we assume that during this part of the evolution mass 
and orbital angular momentum are conserved, the only free pa- 
rameter is the initial mass ratio qu = Mn/Mn. For each of the 
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Fig. 12. Results of the spiral-in calculations with period limits 
for conservative mass transfer as in Fig.[^ but for all systems. 
The number in the upper left corner of each panel is the number 
of systems that lie between the period limits. 
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Fig. 13. Results of the spiral-in calculations (grey symbols), 
obtained as in Fig.^2 an d the solutions of calculations of 
conservative evolution (black symbols). Only the six systems 
shown have spiral-in solutions within the period limits (see 
Fig.1121. The numbers in the lower left and lower right corners 
are the numbers of plotted spiral-in solutions and conservative 
solutions respectively. 

spiral-in solutions we selected, we chose five different values 
for ^n, evenly spread in the logarithm: 1.1, 1.3, 1.7, 2.0 and2.5. 
The total number of conservative models that we calculated is 
570, of which 270 resulted in a double white dwarf. The ma- 
jority of the rest either experienced dynamical mass transfer or 
evolved into a contact system. A few models were discarded 
because of numerical problems. The results of the calculations 
for the conservative first mass transfer are compared to the so- 
lutions of the spiral-in calculations in Fig. 1131 

The systems that result from our conservative models gen- 
erally have longer orbital periods than the intermediate systems 
that we are looking for. This means that stable mass transfer 
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in the models continues beyond the point where the desired 
masses and orbital period are reached. The result is that M\ m 
is too small and that M2 m and P m are too large. The reason 
that mass transfer continues is that the donor star is not yet 
sufficiently evolved: the helium core is still small and there is 
sufficient envelope mass to keep the Roche lobe filled. White 
dwarfs of higher mass would result from larger values of q\\. 
This way, the initial primary is more massive and the initial pe- 
riod is longer, so that the star fills its Roche lobe at a slightly 
later stage in evolution. Both effects increase the mass of the 
resulting white dwarf. However, if one chooses the initial mass 
ratio too high, the system evolves into a contact binary or, for 
even higher qy u mass transfer becomes dynamically unstable. 
In both cases the required intermediate system will not be pro- 
duced. These effects put an upper limit to the initial mass ratio 
for which mass transfer is still stable, and thus an upper limit 
to the white-dwarf mass that can be produced with stable mass 
transfer for a given system mass. Our calculations show that 
conservative models with an initial mass ratio of 2.5 produce 
no double white dwarfs. Apparently this value of q\\ is beyond 
the upper limit. The solutions in Fig.^^with a final mass ratio 
close to or in agreement with the observations come predomi- 
nantly from the models with initial mass ratios of 1.7 and 2.0. 

Because small deviations in the masses and orbital period 
of the intermediate systems can still lead to acceptable double 
white dwarfs, we monitor the evolution of these systems to the 
point where the secondary fills its Roche lobe and determine 
the mass of the second white dwarf Mn from the helium-core 
mass of the secondary at that point. Because the secondary in 
the intermediate binary is slightly too massive in most cases, it 
is smaller at a given core mass (see Fig.^ so that the mass of 
the second white dwarf becomes larger than desired. Combined 
with an undermassive first white dwarf this results in a too large 
mass ratio qn . This is shown in Fig.^J where the values for q^f 
for our conservative models are compared to the observations. 
The Figure also shows the difference in age of the system be- 
tween the moment where the second white dwarf was formed 
and the moment when the first white dwarf was formed (At). 
This difference should be similar to the observed difference in 
cooling age between the two components of the binary (see 
Table^. The vertical dotted lines show this observed cooling- 
age difference with an uncertainty of 50%. 

FigurefT^lshows that of the six systems presented, only two 
have a mass ratio within the observed range, although values 
for the other systems may be close. We see that the mass ra- 
tios of the solutions for most of the systems are divided in two 
groups and the difference in mass ratio can amount to a factor 
of 2 between them. The division arises because in most models 
the common envelope is supposed to occur on the short giant 
branch of stars that are more massive than 2M & . If the sec- 
ondary is slightly smaller and the orbital period slightly longer 
than it should be, the star can ignite helium in its core and start 
shrinking before it has expanded sufficiently to fill its Roche 
lobe. When this star expands again after core helium exhaus- 
tion, it has a much more massive helium core and produces a 
much more massive white dwarf than desired (see Fig.^. Thus, 
a small offset in the parameters of the model after the first mass- 
transfer phase can result in large differences after the spiral-in. 
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Fig. 14. The mass ratio of model double white dwarfs formed 
by a conservative first mass transfer and a common envelope 
with spiral-in, against the age difference between the two com- 
ponents. The dashed horizontal lines show the observed range 
of possible mass ratios as shown in Tabled The dotted verti- 
cal lines are the estimated cooling-age differences +50% (see 
Tabledi. 

Of the 270 stable models shown in Fig.[H| 126 (47%) are in 
the group with lower mass ratios (<?2f ~ 1-7). 

The modelled mass ratios for the systems WD 0957-666 
and WD 1101+364 are close to the observed values, and we 
find that this is especially true for the models on the low-mass 
end of the range in observed white-dwarf masses we used. 
This can be understood, because the maximum mass of a white 
dwarf that can be created with conservative mass transfer is 
set by the total mass in the system. The system mass is deter- 
mined by the spiral-in calculations in Sect. 15.21 where we find 
that the total mass that is available to create these two systems 
lies between about 2 and 3.5 M . This system mass is simply 
insufficient to create white dwarfs with the observed masses. 
If we would extend the uncertainty in the observed masses 
to allow lower white-dwarf masses, it seems likely that we 
could explain these two double white dwarfs with a conserva- 
tive mass-transfer phase followed by a common envelope with 
spiral-in. The same could possibly be achieved with stable, 
non-conservative mass transfer. If the mass were transferred to 
the accretor and subsequently partially lost from the system in 
an isotropic wind, this would stabilise the mass transfer. The 
mass transfer would then still be stable for slightly longer ini- 
tial periods, so that higher initial primary masses are allowed. 
Both effects result in higher white-dwarf masses. 

All 126 stable solutions in the lower group of mass ratios 
(<?2f ~ 1-7) have a ce > 1 and 83 (66%) have c ce < 5. If we 
become more demanding and insist that a ce should be less than 
2, we are left with 14 solutions, all for WD 0957-666. These 
solutions all have a ce > 1.6. If we additionally demand that 
the age difference of these models be less than 50% from the 
observed cooling-age difference, only 6 solutions are left with 
age differences roughly between 190 and 410 Myr, a ce: > 1.8 
and 1.32 < q 2 t < 1.44. 

We conclude that although the evolutionary channel of con- 
servative mass transfer followed by a spiral-in can explain 
some of the observed systems, evolution along this channel 
cannot produce all observed double white dwarfs. We must 
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therefore reject this formation channel as the single mechanism 
to create the white-dwarf binaries. The reason that it fails to 
explain some of the observed white dwarfs is that the observed 
masses for the first white dwarfs in these systems are too high 
to be explained by conservative mass transfer in a binary with 
the total mass that is set by the spiral-in calculations. Allowing 
for mass loss from the system during mass transfer could re- 
sult in better matches. However it is clear from Fig.^]that this 
will certainly not work for at least 5 of the 10 observed systems 
because their orbital periods are too large. We will need to con- 
sider other prescriptions in addition to stable mass transfer to 
produce the observed white-dwarf primaries for these systems. 

6.2. Unstable mass transfer 

In this section we try to explain the formation of the first white 
dwarf in the intermediate systems shown in Fig.^|by unstable 
mass transfer. Mass transfer occurs on the dynamical timescale 
if the donor is evolved and has a deep convective envelope. 
There are two prescriptions that predict the change in orbital 
period in such an event. The first is a classical common en- 
velope with a spiral-in, based on energy conservation as we 
have used in Sect.|5] T he second pres cription was introduced 
in this context bvlNelemans et alJ J2000T) and further explored 
bv lNelemans & Tout! (12005), and uses angular- momentum bal- 
ance to calculate the change in orbital period. Where the first 
prescription results in a strong orbital shrinkage (spiral-in) for 
all systems, in the second prescription this is not necessarily 
the case, so that the orbital period may or may not change dras- 
tically using the same efficiency parameter, while the envelope 
of the donor star is lost. 

In both scenarios we are looking for an initial binary of 
which the components have masses My and M-^. The primary 
will evolve fastest, fill its Roche lobe and eject its envelope due 
to dynamically unstable mass loss, so that its core becomes ex- 
posed and forms a white dwarf with mass M\ m . We assume that 
the mass of the secondary star does not change during this pro- 
cess, so that M2i = M2m- We use the model stars from Sect. 14.11 
as the possible progenitors for the first white dwarf. The orbital 
period before the envelope ejection is again determined by set- 
ting the radius of the model star equal to the Roche-lobe radius 
and applying Eq.|5] where the subscripts 'm' must be replaced 
by 'i'. 

Because we demand that M\\ > Mi\, the original secondary 
can be any but the most massive star from our grid and the total 
number of possible binaries in our grid is Y}„=\ n ~ 19701 for 
each system we want to model. The total number of systems 
that we try to model is 121: the 11 observed systems (the 10 
from Tablefflplus the system WD 1704+481b) times 1 1 differ- 
ent assumptions for the masses of the observed stars (between 
+0.05 M from the observed value). We have thus tried slightly 
less than 2.4 million initial binaries to find acceptable progen- 
itors to these systems. All these possible progenitor systems 
have been filtered by the following criteria, in addition to the 
ones already mentioned in Sect.|5] 

1 . the radius of the star is larger than the radius at the base of 
the giant branch R > Rbgb, which point is defined by Eq.Q 



2. the mass ratio is larger than the critical mass ratio for dy- 
namical mass trans fer q > q cl -i t as defined by Eq. 57 of 
iHurlev et alJ §002). Together with the previous criterion, 
this ensures that the mass transfer can be considered to pro- 
ceed on the dynamical timescale, 

3. the time since the ZAMS after which the first white dwarf is 
created ti is less than the same for the second white dwarf 
(T2) and, additionally, T2 < 13Gyr. 

After we filter the approximately 2.4 million possible pro- 
genitor systems with the criteria above, about 204,000 sys- 
tems are left in the sample (8.5%) for which two subsequent 
envelope-ejection scenarios could result in the desired masses, 
provided that we can somehow explain the change in orbital 
period that is needed to obtain the observed periods. For each 
of the two prescriptions for dynamical mass loss we will see 
whether this sample contains physically acceptable solutions 
in the sections that follow. 

6.2.1. Classical common envelope with spiral-in 

The treatment of a classical common envelope with spiral-in 
based on energy conservation has been described in detail in 
Sect. [5] and therefore need not be reiterated here. In the cal- 
culations described above, Eq.@] provides us with the param- 
eter a ce i for the first spiral-in. In order to use Eq.0]the sub- 
scripts 'm' must be replaced by 'i' and the subscripts 'f ' by 
'm'. The values of the common-envelope parameter for the 
first spiral-in must be physically acceptable and we demand 
that 0.1 < Ccei < 10. When we apply this criterion to the re- 
sults of our calculations, only 25 possible progenitors out of 
the 204,000 binaries in our sample survive. All 25 survivors 
are solutions for WD 0135-052 and have a ce i ^ 2.5. 

We find that of the systems that pass the criterion in the sec- 
ond spiral-in and have 0.1 <a C e2 ^ 10, most (99%) need a nega- 
tive o-cei in order to satisfy Eq.0] so that we reject them. We can 
clearly conclude that the scenario of two subsequent classical 
common envelopes with spiral-in can be rejected as the for- 
mation scenario for any of the ob served double white dw arfs. 
This confirms the conclu sions of Nelema ns et alJ (|2000) and 
iNelemans & Toutl J2005h . based on the value of the product 
a ce /i env , where /l env is the envelope-structure parameter defined 
in Eq.|2l 

6.2.2. Envelope ejection with angular-momentum 
balance 

The idea to determine the change in orbital period in a common 
envelope from balance o f angu lar momentum originates from 
[ Paczvnski & Ziolkowskild 1967b . In lNelemans et alJ J200(A and 
Nelemans & Tout (2005) the prescription was used to model 
observed double white dwarfs. The principle is similar to that 
of a classical common envelope, here with an efficiency pa- 
rameter that we will call y in the general case. In this section 
we will use three slightly different prescriptions for mass loss 
with angular-momentum balance requiring three different def- 
initions of y. For all three prescriptions the mass loss of the 
donor is dynamically unstable and its envelope is ejected from 



M.V. van der Sluys et al.: Modelling the formation of double white dwarfs 



15 



the system. Because not all of these prescriptions necessarily 
involve an envelope that engulfs both stars, we shall refer to 
them as envelope ejection or dynamical mass loss rather than 
common- envelope evolut i on. Th e first prescription is that de- 
fined by iNelem ans et al. (2000), where a common envelope 
is established first, after which the mass is lost from its sur- 
face. The mass thus carries the average angular momentum of 
the system and we will call the parameter for this prescription 
y s . In the second prescription the mass is first transferred and 
then re-emitted with the specific orbital angular momentum of 
the accretor. We will designate y a for this prescription. In the 
third prescription the mass is lost directly from the donor in 
an isotropic wind and the corresponding parameter is ya- We 
will call the companion to the donor star 'accretor', even if no 
matter is actually accreted. 

The prescription of dynamical mass loss with the average 
specific angular momentum of the initial system for the first 
mass-transfer phase, using this and earlier subscript conven- 
tions, is: 



Prescription j\ 



7o,i Jilocti 



Ji 



r.si 



My, - M lm 



(8) 



where J is the total orbital angular momentum ( INelemans et alJ 
2000). Our demands for a physically acceptable solution to ex- 
plain the observed binaries is now 0.1 < y s i < 10 for the first 
envelope ejection and 0.1 < a ce2 < 10 for the second. From the 
set of about 204,000 solutions we found above, almost 150,000 
(72%) meet these demands and nearly 134,000 solutions (66%) 
have values for y s i between 0.5 and 2, in which all observed 
systems are represented. 

We tried to constrain the ranges for y s i and a ce2 as much 
as possible, thereby keeping at least one solution for each ob- 
served system. We can write these ranges as (yo - To + ^) 
and (ao - ^,ao + 4p), where yo and ao are the central val- 
ues and Ay and Ao- are the widths of each range. We indepen- 
dently varied yo and ao and for each pair we took the smallest 
values of Ay and Aa for which there is at least one solution 
for each observed system that lies within both ranges. The set 
of smallest ranges thus obtained is considered to be the best 
range for y s i and a ce2 that can explain all systems. Because it 
is harder to trifle with the angular- momentum budget than with 
that of energy, we kept the relative width of the range for y s i 
twice as small as that for a ce2 (2— = —). Our calculations 
show that changing this factor merely redistributes the widths 
over the two ranges without affecting the central values much 
and thus precisely which factor we use seems to be unimpor- 
tant for the result. We find that the set of narrowest ranges that 
contain a solution for each system is 1.45 < y s i < 1.58 and 
0.61 <a C e2<0.72. These results are plotted in Fig. 1151 

We can alternatively treat the second envelope ejection with 
the angular-momentum prescription as well, where we need 
to introduce a factor y S 2 by replacing all subscripts 'm' by 
'f and all subscripts 'i' by 'm' in Eq.[8] Again we search 
for the narrowest ranges of y s i and y S 2 that contain at least 
one solution per observed system. We now force the relative 
widths of the two ranges to be equal. The best solution is then 
1.16<y s i< 1.22 and 1.62 <y s2 < 1.69. 
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Table 2. Narrowest ranges for y and a ce that contain at least one 
solution to the envelope-ejection scenario per observed system 
and their central values. The six different prescriptions are ex- 
plained in the main text. 



In both prescriptions above (y s io , C e2 and y s iy S 2) we find 
that the values for y lie signific antly above unity . This is in 
accordance with t he fin dings of INelemans et all J2000) and 
Nelemans & Tout (2005), but slightly discomforting because 
there is no physical reason as to why y should have this, or in- 
deed any other, particular value. We will therefore rewrite Eq.|S] 
for the case where the mass is lost with the specific orbital an- 
gular momentum of one of the stars in the binary, so that we 
can expect that y » 1 . In order to do this we use the equations 
derived by Sober man et alJ \ 19971) in their Section 2.1. They 
define the fractions of mass lost by the donor: a w the fraction 
that is lost in an isotropic wind from the donor, yS w the fraction 
that is transferred and then lost in a wind by the companion 
and e w the fraction that is accreted by the companion, where 
we introduced the subscript 'w' to avoid confusion with a ce . 
We assume that no matter is accreted, so that a w + /? w = 1 and 
e w = and ignore the finite size of the star by putting A w = 1 . 
Their Eq. 24 then gives (replacing their notation by ours): 

"*w -I 

1 +qi 



J 111 
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exp [/?„ (q m - qi)] . 



(9) 



where we will consider the cases where a w =0; /3 V — 1, describ- 
ing isotropic re-emission by the accretor, and ar w = 1 ; yS w = 
for an isotropic wind from the donor. Their q is defined as 
'Wdonor/'Maccretor- We can now replace Eq.|H]with similar expres- 
sions for these two cases: 
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(11) 



where we introduced the parameters y a i and yai to describe 
deviations from the idealised case of purely isotropic mass loss. 
We treat y a i and y<ji as free parameters, but expect that y « 1 in 
the picture that mass is preferentially lost from either donor or 
accretor. We used M tot = M\ +M 2 to denote the total mass of the 
system. The results of the analysis described above, but now for 
the modified definitions of y, for the ya and yy scenarios, each 
with a w = 0, /? w = 1 (re-emission from the accretor, y a ~ 1) 
and a „ — 1 , /3 W = (wind from the donor, y<j « 1) are shown in 
Tableland compared to the previous results. 

We see that the values for y change drastically, as may be 
expected. The fact that the values for a ce change slightly has 
to do with the fact that we now select different solutions to the 
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Fig. 15. The distribution of the logarithm of the envelope-ejection parameters for solutions of the double dynamical mass-loss 
scenario. Each dot represents a system that evolves through an episode of dynamical mass loss with y s i and then a common 
envelope with spiral-in with a c< a to form one of the observed white dwarfs. Only solutions for which 0.1 < y s i < 10 and 
0.1 < a ce 2 < 10 are shown. Hence, each dot in one panel has a counterpart in the other panel. The dashed vertical lines show 
y s i = 1.45, 1.58 in the left panel and ar ce 2 = 0.61, 0.72 in the right panel. Black dots in one panel indicate solutions that lie between 
the dashed lines in the other panel; black dots between black lines therefore lie in the preferred range in both panels. The vertical 
position of each dot within the line for its system shows the deviation from the observed secondary mass Mi. From the lower to 
the upper dots in a line these masses are: Mh = M? - 0.05, M? - 0.04, . . . , Mi + 0.05 M Q . The total number of solutions in the 
Figure is nearly 150,000. The Figure is made after Nelemans & ToutH2005l) . 



calculations than before. The prescription y a iTa2 has very nar- 
row ranges. Numerically, the prescription ydiQ'ce2 seems very 
attractive: ya » 1.0 and a ce 2 ~ 0.6. Although the value for 
a ce 2 is lower than unity, it may not be unrealistic that 40 % of 
the freed orbital energy is emitted by radiation. This is the sce- 
nario where the mass is lost in an isotropic wind by the donor 
in the first dynamical mass-loss episode and the second mass 
loss is a canonical common envelope with spiral-in. We also 
notice that in order to get a sufficiently strong spiral-in from 
an envelope ejection by a strong donor wind (the 7diTd2 pre- 
scription), the mass has to be lost with almost three times the 
specific angular momentum of the donor. 

6.3. Formation by multiple prescriptions 

So far, we assumed that all ten observed double white dwarfs 
were formed by one and the same prescription. Although some 
prescriptions are clearly better in explaining the formation of 
all the observed systems than others, none of them is com- 
pletely satisfactory, mainly because the parameters y or a ce are 
far from the desired values. Furthermore, there is no reason 
why the ten systems should all have evolved according to the 
same prescription in nature if there are several options avail- 
able. We therefore slightly change our strategy here by assum- 
ing that different envelope-ejection prescriptions, described in 
Sect. 16.2.21 can play a role in the formation of each of the ob- 
served systems. 

For the dynamical mass loss, we now demand that y and 
or ce are close to unity. Because angular momentum should be 
better conserved than energy, we accept solutions with 0.90 < 
y a> d < 1.10 and 0.80 < ar ce < 1.20. For the prescription by Eq.[3] 



iNelemans & Toutl J2005I) show that all systems can be ex- 
plained with 1.50 < y s < 1.75, which we adapt to 1.46<y s < 1.79 
to give it the same relative width. The total number of solutions 
between these ranges is about 42500, divided over 6 prescrip- 
tions and 1 1 systems (on average 3.5% of the solutions selected 
in Sect. 16.2^ For each observed system and each prescription, 
we look whether there is at least one solution with a envelope- 
ejection parameter within these ranges. The results are shown 
as the first symbol in each entry of Table[3] The plus signs show 
which prescription can explain the mass ratio of an observed 
double white dwarf with envelope ejection parameters within 
the chosen ranges of y and a ce . The table confirms that the pre- 
scription y s y s can explain all observed systems in this sense. 
The prescription y^y-d can do the same if we are satisfied with 
a solution for WD 1704+481 a or b. We did not include the 
columns for y- d yd and ydTd, because they cannot explain any 
of the observed systems. This can readily be understood, since 
in these two scenarios the last envelope is lost by a wind from 
the donor. This process will usually widen the orbit, rather than 
cause the spiral-in that is needed to explain the observed bina- 
ries. We have to expand the y-ranges to 0.25-1.75 in order to 
get the first solution for just a single system with one of these 
two prescriptions and only by allowing values for y up to 2.8, 
we can find a solution for each binary. 

6.4. Constraining the age difference 

The large number of solutions found in the previous section al- 
lows us to increase the number of selection criteria that we use 
to qualify a solution as physically acceptable. We now include 
the age difference of the components in our model systems and 
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System 


1 : 7s«ce 


2: 7s7s 


3: yaQ-cc 


4: 7a7a 


5: 7d»cc 


6: 7d7a 


Opt. res. 


Best prescr. 


0135 


-/- 


+/~ 


+/~ 


-/- 


+/~ 


+/~ 


+/~ 


2,3,5,6 


0136 


+/+ 


+/+ 


+/~ 


+/~ 


+ /+ 


+/+ 


+ /+ 


1,2,5,6 


0957 


+/+ 


+/+ 


-/- 


+/- 


+ /+ 


+/+ 


+/+ 


1,2,5,6 


1101 


+/~ 


+/- 


+/- 


-/- 


+/~ 


+/~ 


+/~ 


1,5,6 


1115 


+/~ 


+/+ 


+/~ 


+/~ 


+ /+ 


+/+ 


+ /+ 


2,5,6 


1204 


-/- 


+/- 


+/- 


+/- 


+/- 


+/+ 


+ /+ 


6 


1349 


+/+ 


+/+ 


+/+ 


+/+ 


+/+ 


+/+ 


+/+ 


1,2,3,4,5,6 


1414 


-/- 


+/+ 


-/- 


+/+ 


-/- 


+/+ 


+/ + 


2,4,6 


1704a 


+/- 


+/- 


-/- 


-/- 


-/- 


-/- 


+/- 


1,2 


1704b 


+/- 


+/- 


-/- 


+/- 


+/- 


+/- 


+/- 


1,2,4,5,6 


2209 


+/+ 


+/+ 


-/- 


-/- 


+/~ 


+/+ 


+/+ 


1,2,6 



Table 3. Comparison of the different prescriptions used to reconstruct the observed double white dwarfs. The symbols +, ~ and 
- mean that the model solutions are in good, moderate or bad agreement with the observations. The first of the two symbols in 
each column is based on the mass ratio only and the second includes the age difference. The method for obtaining the first symbol 
in each entry is described in Sect. 16.31 that for the second symbol in Sect. 16.41 The symbols in the headers of the columns labelled 
1-6 are explained in the main text. The columns for y a yd and yajd were left out because they do not contain any solutions. The 
last two columns show the optimum result and the prescriptions that give this result (1-6). For the different prescriptions we 
demanded that 1 .46 < y s < 1 .79, 0.90 < y a < 1 . 10, 0.90 < y d < 1 ■ 10 and 0.80 < a ce < 1 .20. 



demand that it is comparable to the observed cooling-age dif- 
ference for that system. The age difference in our models is 
the difference in age at which each of the components fills its 
Roche lobe and causes dynamical mass loss. In addition, we 
reject the solutions for which the initial mass ratio qn < 1.03. 
In these binaries, both components have probably evolved be- 
yond the main sequence and it is uncertain what would happen 
to the secondary at the envelope ejection of the primary. 

Tabled lists the number of model solutions for each pre- 
scription and each system. The columns labelled 1-6 represent 
the same prescriptions as those columns in Table[3] The first 
number in each of these columns is the number of solutions 
that is found within the same ranges for y and a ce as we used in 
Table|5] This means that a minus sign in that table corresponds 
to a zero in Tabled Behind the entries with a positive num- 
ber of solutions the range of age difference that these solutions 
span is shown. Again, the columns for y- d yd and ydTd are not 
displayed, because they do not contain any solutions for any 
system. All prescriptions that are listed in Table|4]provide a so- 
lution for more than one observed system, and each observed 
system has at least one prescription that provides it with a solu- 
tion. The number of solutions per combination of prescription 
and observed system ranges from zero to several hundreds and 
the age differences of the accepted models lie between 36 Myr 
and more than 12 Gyr. 

We will use Table|4] to compare the age differences of the 
models to the observed values and use this comparison to judge 
the 'quality' of the model solutions. We will assume that if 
the age difference in the model lies within 50% of the mea- 
sured cooling-age difference (the range in the second column 
of Table|4} that this is a good agreement which we will assign 
the symbol '+'. If the difference is larger than that, but smaller 
than a factor of five we will call it 'close' and assign a 
Cases where the nearest solution has an age difference that is 
more than a factor of five from the observed value is considered 
'bad' and assigned the symbol '-'. If we do this for all cases, 



we obtain the second symbol for each entry in Tabled which 
we can use to directly compare the quality of the solutions for 
each prescription and each observed system. 

If we compare the first and second symbol of each entry in 
Table|5J we see that imposing the extra demand of a model age 
difference (At) that is comparable to the observations rejects a 
significant number of solutions. We find that from the 42500 
solutions that we found for the six prescriptions, based on final 
mass ratio only, 39400 are left if we demand that q\\ > 1.03, 
3 1 100 if we additionally demand that the model value of At lies 
within a factor of 5 from the observed value, and 26800 if we 
demand that At given by the model lies less that 50% from ob- 
served. This means that the plusses in some entries in the Table 
now become tildes or minuses, because no model solution can 
be found with an acceptable age difference for the given combi- 
nation of observed system and envelope-ejection prescription. 
This is especially the case for WD 0135-052, WD 1101+364 
and WD 1704+481 and to a lesser extent WD 1204+450. 

The prescription ydy a now gives slightly better results than 
y s y s , and we find that ydy a gives the same or a better result 
than y s ys for each system except WD 1 704+48 1 a, for which the 
prescription ydy a does not provide a solution with the correct 
mass ratio. The envelope-ejection prescription ydy a can explain 
all systems except WD 1704+481 fairly well, as we shall see 
in the next section that the two tildes for WD 0135-052 and 
especially WD 1101+364 are very close to the demands for a 
plus. FigurefTBI shows the solutions for the yd 7a scenario for 
each system schematically. 

6.5. Description of selected solutions 

The goal of this research is, of course, to find out whether we 
can somehow explain the formation of the observed double 
white dwarfs. If this is the case, we hope to learn firstly which 
prescriptions govern this formation and secondly what the pro- 
genitor systems are that evolve to the observed white-dwarf bi- 
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System Obs. At Number of solutions and model age differences (Myr) 





(Myr) 


1: ysQ-ce 


2: TsTs 


3: y a a C e 


4: y a 7a 


5: rdttce 


6: TdTa 


0135 


175-525 





5, 636-1138 


44, 1454-3598 





775, 1198-3608 


57, 53-1495 




??S 67S 


jsn 170 IIS 


407 170-Q71 


?fll 1^(5fi 7^08 


7 76 1^5(W10(S0 


7?7? 170 540^ 


700? 161 ^8^8 

j 7(7 J , 1U1 — JOJO 


0957 


1 uj too 


1R4 764 107Q 

JOT, — 1U / 7 


547? 1 5071 


Q 


5 ?04fi-^(S50 

-7, Z.7 t tU JUJ7 


774 17(S 1070 


755? 00 (S??(S 

7-7-7Z, 7U U^XU 


1 101 


108 ^7^ 


747? 8(S8 1 1^31 

iTii, OUO — ±1.7.71 


4070 17S^ 11(S14 


0000-17755 





7 ? 000 1 087 
j j, yy\j 1 ;70 / 


74 ^fi^ ?6^(S 


1115 


80 740 


6?5 744. 1141 


?07? 70S 1075 


70 648 1 170 


4 S7^ 50<S 


740-1^44 


99Q ^09 

J OL/, 7 — J7Z 


1204 


40-120 




5 808-1087 




4 1 307-1 505 


7 57 1716-5670 


42, 74-1495 


1349 




767, 215-760 


633, 170-840 


S57, 1690-11521 


274, 1356-6374 


7892, 215-6505 


4283, 157-6168 


1414 


100-300 





74, 83-146 





49, 290-920 





6, 167-207 


1704a 


-30- -10 


27, 858-1294 


770, 52-2046 














1704b 


10-30 


6, 465-553 


7 6 3, 307-933 





9, 1060-1313 


6, 465-553 


377, 205-1498 


2209 


250-750 


6, 552-967 


7049, 61-967 








2, 933-975 


767, 61-1060 



Table 4. Results for the various evolution scenarios for double white dwarfs with two unstable mass-transfer episodes. The range 
of observed At is the observed cooling-age difference + 50%. Columns labelled 1 through 6 give the number of model solutions 
for each scenario followed by the range in age difference of these solutions in Megayears. The columns with y^ya and ydjd were 
left out, because they do not contain any solutions. In addition to a value for y or a in the desired ranges (see Tabled, we demand 
that q n > 1.03. 



0135 

0136 

0957 

I 10: 

1115 

1204 

1349 

1414 

1704a 

1704b 

2209 



-1 



-0.5 







0.5 



0135 
0136 



1414 
1704a 
1704b 
2209 



-1 



-0.5 





l0 § ? a 3 



0.5 



Fig. 16. Solutions for the y^Ta scenario for each system. The representation is the same as in Fig.[21 The vertical dashed lines 
are at y = 0.9 and y — 1.1 in both panels. Black dots are solutions that lie between these lines in the other panel. Black dots that 
lie between the lines thus have 0.9 < yai ^ 1-1 and 0.9 < y a 2 < 1.1. The total number of solutions exceeds 200,000. 



naries. Although we do not find one prescription that can ex- 
plain all observed systems in a perfectly satisfying way, the pre- 
scription yaTa comes very close. In this section we present the 
evolution of some of the best solutions among our calculations. 
We list the five main parameters that describe the evolutionary 
scenario of a solution (two initial masses, the initial period and 
the two envelope-ejection parameters) and two minor ones (the 
actual final masses that the models give) in a table. Because we 
present solutions for six different formation prescriptions and 
eleven observed systems that may have more than one 'best' 
solution this table is quite large. This is particularly the case 
because we want to remove the arbitrarily chosen boundaries 
that we have used so far to qualify a solution. We therefore list 
at least one solution per prescription per system, independent of 
how far its parameters lie from the preferred values. We chose 
to publish the complete table in electronic form only and give 



an excerpt of it in this paper in Table[5] In this way, the reader 
may verify how particular models do or do not work. 

We manually picked the 'best' solutions for a given com- 
bination of formation prescription and observed system, in the 
sense that the solution has a y close to unity (or, in case of 
y s , close to 1.63), an a ce close to the range of 0.5-1.0 and an 
age difference that is close to the observed value. In the cases 
where there are different solutions that each excel in a differ- 
ent one of these three properties, we may present more than 
one solution. If there are several solutions that are similar on 
these grounds, we prefer those with lower initial masses and 
with qu > 1.03. We then leave it to the reader to judge whether 
these solutions are acceptable. The values for <?2f and Pf are 
identical to the value listed in Tableland we do not list q2t in 
Table|5] The intermediate masses are also left out of the table, 
because no matter is accreted during the dynamical mass loss 
and thus M\ m = Mu and Mi m = Mj\ in our models. The num- 
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Nr. 


WD 


Presc. 


7i 


72, 


At (Myr) 


A(At) 


1"2 


M n 


M 2i 


qn 


Pi 


?2m 


Pm 


Mi, 


M 2 f 


Pf 










ffce2 


Obs 


Mdl 


% 


Gyr 


M e 


M 




d 




d 


Me 


Me 


d 


4 


0135 


Ts7s 


1.75 


1.57 


350 


636 


82. 


1.36 


2.59 


2.06 


1.26 


150.8 


3.56 


110.9 


0.58 


0.52 


1.56 


10 


0135 


7d7a 


1.11 


0.94 


350 


118 


66. 


0.41 


3.30 


2.90 


1.14 


36.28 


6.21 


41.10 


0.47 


0.42 


1.56 


12 


0136 


7s«ce 


1.48 


1.01 


450 


449 


0.2 


3.73 


1.44 


1.38 


1.04 


72.27 


4.16 


265.9 


0.33 


0.42 


1.41 


21 


0136 


7d7a 


0.96 


1.05 


450 


450 


0.0 


2.39 


1.70 


1.59 


1.07 


106.1 


4.36 


371.4 


0.37 


0.46 


1.41 


24 


0957 




1.62 


1.00 


325 


427 


31. 


1.16 


2.34 


2.00 


1.17 


8.110 


6.66 


28.52 


0.30 


0.34 


0.06 


32 


0957 


7d7a 


1.00 


1.01 


325 


317 


2.5 


1.52 


1.98 


1.83 


1.08 


26.17 


5.61 


79.26 


0.33 


0.37 


0.06 


44 


1101 


7d7a 


1.10 


0.98 


215 


322 


50. 


0.74 


2.87 


2.34 


1.23 


22.02 


5.97 


28.23 


0.39 


0.34 


0.14 


46 


1115 


7s«ce 


1.79 


1.00 


160 


239 


49. 


0.50 


3.70 


2.94 


1.26 


1693. 


3.58 


980.4 


0.82 


0.69 


30.09 


57 


1115 


7d7a 


0.97 


1.04 


160 


240 


50. 


0.32 


5.42 


3.42 


1.58 


201.2 


3.84 


1012. 


0.89 


0.75 


30.09 


60 


1204 


7s7s 


1.89 


1.25 


80 


80 


0.5 


0.31 


3.60 


3.21 


1.12 


69.10 


6.07 


39.45 


0.53 


0.46 


1.60 


68 


1204 


7d7a 


1.09 


0.92 


80 


100 


25. 


0.38 


3.34 


2.98 


1.12 


15.47 


6.32 


19.99 


0.47 


0.41 


1.60 


70 


1349 


7s«ce 


1.45 


1.01 


o 


461 


0.0 


4.52 


1.35 


1.32 


1.03 


105.9 


3.77 


364.5 


0.35 


0.44 


2.21 


77 


1349 


7d7a 


0.95 


0.98 


o 


101 


0.0 


1.58 


1.86 


1.81 


1.03 


63.44 


5.19 


241.2 


0.35 


0.44 


2.21 


80 


1414 


7s7s 


1.49 


1.76 


200 


1 12 


44. 


0.87 


2.55 


2.43 


1.05 


525.4 


4.09 


1783. 


0.59 


0.76 


0.52 


85 


1414 


7d7a 


0.95 


0.99 


200 


188 


5.9 


0.43 


3.51 


3.09 


1.14 


70.81 


5.99 


358.3 


0.52 


0.66 


0.52 


89 


1704a 


7s<*ce 


2.05 


0.43 


-20 


7 


135. 


1.36 


2.03 


1.90 


1.07 


252.8 


3.51 


96.02 


0.54 


0.38 


0.14 


98 


1704a 


7d7a 


1.11 


1.13 


-20 


52 


360. 


1.41 


2.06 


1.88 


1.09 


40.37 


3.66 


65.66 


0.51 


0.36 


0.14 


101 


1704b 


7s«ce 


1.74 


0.76 


20 


285 


1326. 


0.75 


2.76 


2.55 


1.08 


49.12 


6.40 


107.3 


0.40 


0.57 


0.14 


107 


1704b 


7d«ce 


1.03 


0.15 


20 


182 


810. 


2.23 


1.68 


1.65 


1.01 


212.1 


4.08 


478.6 


0.41 


0.58 


0.14 


113 


2209 


7s7» 


1.62 


1.69 


500 


477 


4.5 


1.37 


2.40 


2.00 


1.20 


165.4 


3.64 


258.7 


0.55 


0.55 


0.28 


122 


2209 


7d7a 


1.04 


1.05 


500 


340 


32. 


0.50 


4.15 


2.94 


1.41 


98.45 


4.67 


294.3 


0.63 


0.63 


0.28 



Table 5. Selected model solutions for the double envelope-ejection scenario. This table is an excerpt of the total list of 122 entries. 
The first nine columns show the number of the entry, the double white dwarf that the model is a solution for, the prescription used, 
the two envelope-ejection parameters, the age difference of the two components as observed and in the model (At) in Myr, the 



relative difference between the observed and model age difference, defined as A(At) = | t "'a t b r " bs | m %> me tmie °f me formation 
of the double white dwarf since the ZAMS (j2) in Gyr. The last nine columns list binary parameters: the initial (ZAMS) masses, 
mass ratio and orbital period, the intermediate mass ratio and period and the final masses and period. 



bering of the solutions in the excerpted table is the same as in 
the complete version. 

The complete table contains 122 solutions. The initial bi- 
naries have primary masses between 1.27 M@ and 5.42 M . 
Of the 122 solutions, 49% have an initial primary mass less 
than 2.5 M Q and 93% of the primaries are less massive than 
4M Q . Thus, the models suggest that the double white dwarfs 
are formed by low-mass stars, as may be required to explain 
the observed numbers of these binaries. Of the initial systems, 
89% have orbital periods between 10 and 1000 days. All pro- 
posed solutions undergo a first envelope ejection described by 
angular-momentum balance of some sort, which allows the or- 
bital period to increase during such a mass-transfer phase. In 
80% of the selected solutions this is the case, and for 61% of 
the solutions the intermediate orbital period is twice or more as 
long as the initial period. Of the 122 solutions, 57% have initial 
mass ratios qu < 1.15 while for 17% qn > 1.5. We included 
only 3 solutions with qn < 1.03. It could be that these initial 
systems evolve into a double common envelope, where the two 
white dwarfs are formed simultaneously and the second white 
dwarf is undermassive. On the other hand, because the orbital 
period increases in most of the first envelope ejections, the out- 
come of such a common envelope is uncertain. We should treat 
these solutions with some scepticism. 

We now briefly discuss the individual solutions that are 
listed in the excerpted Table|5] For WD 0135-052 we find so- 



lutions with acceptable y's or with acceptable At, but not with 
both. In solution 4 the age difference differs 82% from observed 
while the values fory s lie less than 10% from 1.63. Solution 10 
(yd/a) has a value for yd just outside the 10% region (10.8%) 
and a value for At that lies just outside the 50% region (66%). 
These two solutions therefore have properties that lie not too far 
from the preferred values. Solutions 12 (y s a ce ) and 21 (ydy a ) 
for WD 0136+768 have envelope-ejection parameters within 
10% and 5% respectively from the preferred values, while the 
model ages hardly differ at all from the observed values. Both 
solutions are also a solution for y s y s (in the complete table). 
In solutions 24 and 32 for WD 0957-666 we find envelope- 
ejection parameters with less than 1 % from the preferred val- 
ues and the age differences well within 50% from observed. 
For WD 1101+364, solution is very close to the requirements 
for a plus sign in Table|3] At lies just less than 50% from ob- 
served, however yd = 1.103, so that it falls just outside the 
10% range from the preferred value. Solutions 46 and 57 for 
PG 1115+116 have envelope-ejection parameters 10% and 5% 
respectively from the desired values and acceptable values for 
At. The initial masses are high for these solutions, in accor- 
dance with the fact that these stars are required to form white 
dwarfs with masses as high as 0.8 M . 

WD 1204+450 has few solutions that have an age differ- 
ence as small as observed. Solutions that do have the observed 
age difference, as solution 60, often have values for y or a 
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that lie far from desired. The prescription y&y^ is the only one 
that provides solutions that are qualified as 'good' in Table|5] 
Solution 68 has values for y within 10% from unity and an age 
difference 25% from the observed value. For WD 1349+144 
the cooling ages are not known, although the similar Balmer 
spectra of the two components jKarl et all2 003a) seem to sug- 
gest that At is small. Solutions 70 and 77 give solutions with 
envelope-ejection parameters slightly more that 10% and 5% 
from the preferred values respectively and have age differences 
of 461 and 101 Myr. The complete table shows that in addition 
to this solutions with values for At such as 129, 215, 384, 975 
and 1776 Myr can be found, so that they span a large range 
within which the actual age difference is likely to lie. The pre- 
scriptions y s y s and y^Ta gi ve good solutions for the system 
HE 1414-0848. The solutions 80 and 85 have values for y that 
lie within 10% and 5% respectively from the desired values, 
while their age differences are 44% and 6% from observed. 

Since the observed age difference of WD 1704+48 la is 
-20 Myr, we introduced a system with the reversed mass ra- 
tio (WD 1704+481b) and hence an age difference of +20 Myr. 
Interestingly enough, the solutions with closest age difference 
for WD 1704+48 lb have At> 1 80 Myr, a factor of nine or more 
than observed, as is the case for solution 107, though it has 
a ce = 0.15. Solution 109 has nice values for y, but an age dif- 
ference of 205 Myr. However, for WD 1704+48 la we find solu- 
tions with reasonable envelope-ejection parameters (13% from 
the preferred value) and an age difference of around 50 Myr, 
like solution 98, and with parameters that are more off, but with 
an age difference of only 7 Myr as in solution 89. The system 
WD 1704+48 la seems therefore better explained by our mod- 
els than the system with the reverse mass ratio. Because the 
observed cooling-age difference is only in the order of a few 
per cent of the total age of the system (see Tabled, a change 
of 10% in the determined cooling age of one of the two com- 
ponents is sufficient to alter the age difference from -20 Myr 
to +50 Myr, which would bring it in the range that can be ex- 
plained by our models. For HE 2209-1444, we present solu- 
tions 113 and 122, that have envelope-ejection parameters less 
than 5% from the desired values and age differences that lie 
5% and 32% respectively from the observed cooling-age dif- 
ference. 

Summarising, we find that we can explain the mass ratios 
and orbital periods of the ten observed systems well with the 
scenarios y s y s and ydy a - If we also take into account the ob- 
served age difference between the components, the scenario 
7d7a still works well. However, the allowed ranges for y and 
age difference must be extended a little to include WD 0135- 
052 ( rdl = 1.11, A(At) = 66%) and WD 1101+364 (y d] =1.11), 
while the small age difference observed for WD 1704+48 1 can- 
not be explained. 

7. Discussion 

7.1. Comparison to other work 

In this paper we investigate the f ormation scenarios for double 
white dwarfs first put forward bv lNelemans et alJ J2000h . Their 
paper is based on three double white dwarfs and we expanded 



this to the set of ten that has been observed so far. Rather than 
using analytical approximations, we used a stellar evolution 
code to do most of the calculations. One of the advantages of 
this is that we calculate the binding energy of the donor star at 
the onset of a common envelope, so that we can directly cal- 
culate the common-envelope parameter ar oe without the need of 
the envelope-structure parameter A em , that turns out to be far 
from constant during the evolution of a star (see Fig.|4ji. This 
allows us to demand physically acceptable values for a ce . 

The use of an evolution code instead of analytical expres- 
sions obviously gives more accurate values for instance for the 
core-mass - radius relatio n. Our main conclusion s are neverthe- 
less the same as those of Nelema ns et al.l (120001) . even though 
they are based on a larger sample of observed binaries: sta- 
ble, conservative mass transfer followed by a common enve- 
lope with spiral-in based on energy balance cannot explain the 
formation of the observed systems, and neither can the a Cfl a ce 
scenario of two such spiral-ins. We therefore arrive at the same 
conclusion, that a third mass-transfer prescription is needed 
to explain the first mass-transfer phase of these systems and 
we use their envelope-ejection prescription, based on angular- 
momentum balance (Eq.|HJ. 

Nelemans & Tout (2005) use more advanced fits to stellar 
models, but still need the envelope-structure parameter A env so 
that it is difficult to interpret the values they find for the prod- 
uct ff ce /lenv They use the same ten observed double-lined white 
dwarfs as we do, next to a number of single-lined systems. 
They also conclude that a y-envelope ej ection is needed f or the 
first mass transfer and find, like iNelemans et alJ l|2000), that 
all observed systems can be explained by 1.50 < y s < 1.75, 
for both mass-transfer phases. Alternatively the second mass- 
transfer e pisode can be reconstru cted with < a ce: /l env < 4. 
However, Nel emans & T out (2005) do not discuss the coupling 
of the two solution sets for the two phases, e.g. it is not de- 
scribed how many of the solutions with 1.50 < y s i < 1.75 have 
y S 2 in the same range. We introduced slightly different defi- 
nitions for the y-algorithm in Eqs.^|and^2 so that we can 
demand that y is in the order of unity. We find indeed that we 
can explain the observed masses a nd periods with y„, yd ~ 1 .0. 

We add to the treat ment by Nelem ans et alJ J2000I) and 
INelemans & Tbul {2005) in demanding that, in addition to the 
masses and orbital period, the age difference of our models 
must be comparable to the observed value. It turns out that this 
puts a strong constraint on the selection of model solutions for 
all three definitions of y. However, we can still explain most 
systems, although we need mass loss described by both y a and 
yd to do so. 

The description for dynamical mass loss with the specific 
angular momentum of the donor star (Eall 111 is similar to the 
s cenario of a tidally-enhanced stellar wind (Tout & Eggleton 

El I J 
988 b a). In this scenario the mass loss from a (sub)giant due 

to stellar wind increases up to a factor of 150 with respect to 
Reimers' empirical law jReim ersjl 197 51) when th e star is close 
to filling its Roche lobe. lTout & Eggleton ( 1988b) postulate the 
enhanced wind to explain for instance observed pre-Algol sys- 
tems such as ZHer. In this binary the more evolved star is 10% 
less massive than its main-sequence companion, while only fill- 
ing about half of its Roche lobe. 
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iHanl ill 9981) uses this tidally-enhanced stellar wind in his re- 
search on the formation of double degenerates and concludes, 
among others, that his models that include this enhanced stel- 
lar wind give a better explanation of the observed double- 
degenerate binaries than models that do not include it. The 
enhanced mass loss makes subsequent mass transfer due to 
Roche-lobe overflow dynamically more stable. Envelope ejec- 
tion due to dynamical mass loss is then more often prevented 
and binaries evolve to longer orbital periods before the second 
mass transfer, which is then more likely to produce a CO white 
dwarf. Thus, the enhanced-wind scenario increases the ratio of 
CO-helium double white dwarfs to helium-helium binaries. 

Envelope ejection described by Eq^2 is essentially the 
same as the limiting case in which most or all of the enve- 
lope is lost due to an enhanced wind. lTout & Eggleton ( 1988b) 
show that the tidally-enhanced wind can indeed prevent Roche- 
lobe overflow altogether because the envelope is completely 
lost by the wind and the core becomes exposed. Without an 
enhanced wind, this happens for binaries with an initial mass 
ratio of 2 or less only if they have initial periods of more than 
1000 days. When the tidally-enhanced wind is included, core 
exposure without Roche-lobe overflow occurs for these bina- 
ries with initial periods as short as 10-30 days. 

7.2. Alternative formation scenario for massive white 
dwarfs 

In the present research we have assumed that after envelope 
ejection occurs, the core of the Roche-lobe filling giant be- 
comes a helium or CO white dwarf with no further evolution 
other than cooling. However, helium cores that are more mas- 
sive than 0.33 M are not degenerate and those more massive 
than about 0.5 M will burn most of the helium in their cores 
and produce a CO core. If exposed, they are in effect helium 
stars. For helium stars less massive than about 0.75 M the ra- 
dius hardly changes during the helium (shell) burning, but stars 
more massive than that experience a giant phase. This is shown 
in Fig.ll7h. where the radius of a selection of helium-star mod- 
els is plotted as a function of the CO-core mass. For the more 
massive models in the Figure, the stars expand from the order 
of a few tenths of a solar radius to a few hundred solar radii. 
Thus, helium stars with a core mass Mco ~ 0.7 M may and 
those with Mqo ^0.8 M Q must become giants and could fill 
their Roche lobes as a consequence. The black dots in Fig.ll7h 
indicate the maximum radius for a certain helium-star model 
and if we plot the maximum radii of these and other models 
as a function of the total mass of the star, we obtain Fig. 117b. 
This Figure shows that a helium star more massive than about 
0.83 M must evolve through a giant phase (see e.g. Paczyriski 
1971; Habets 1986). 

There are two double white dwarfs in the observed sam- 
ple with M 2 > 0.6 M Q , PG 1115+116 (both components) and 
HE 1414-0848 (the secondary). The evolutionary scenarios in 
Table|5] suggest that all these stars emerge from the envelope 
ejection with a CO core, except solution57 for PG 1 1 15+116, 
where the 5.42 M primary progenitor produces a 0.89 M he- 
lium core before helium ignites. The Roche-lobe radius of the 
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Fig. 17. Upper panel (a): The radius of a helium star as a func- 
tion of its CO-core mass, for a selection of 15 models with total 
masses between 0.41 and 1.43 M Q . The dots show where the 
maximum radii are obtained and are used for the lower panel. 
The dashed line is the Roche-lobe radius for the intermediate 
primary of PG 1115+116 according to our solution 57. Lower 
panel (b): The maximum radius of a low-mass helium star as 
a function of its total mass, for a selection of 33 models with 
masses between 0.33 and 1.4M . The dots are the data points, 
the solid line connects them to guide the eye. 



0.89 M Q helium star in the intermediate binary is 187 7? ac- 
cording to this solution and shown as the dashed line in Fig.1171 
The Figure shows that the mass and Roche-lobe radius of 
this star are in the proper range to fit the helium-giant scenario. 
We show a simple numerical example to illustrate this scenario. 
The dot in Fig.[T7Ji at M CO =O.88M and R = 171 R Q , just be- 
low the dashed line, is the point where the model of 0.93 M 
from our grid of helium-star models reaches its largest radius. 
The star thus has an envelope mass of only 0.05 M and with 
a ma ss ratio of almost 4, mass transfer would be stable (Eq. 57 
of lHurlev etaD ll2002l) \ If we assume that this star would be 
the primary of solution 57 in Table|5] and that 0.04 M would 
be transferred conservatively, the orbital period after the mass 
transfer would be 1 1 15 d, so that the period would not change 
drastically and the ensuing second envelope ejection would be 
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similar to the one found in solution 57. If the mass were lost 
in a wind, which could be triggered by the fact that the star 
expands, but for which the Roche lobe need not be filled, the 
orbital period would change less than 2% to 1031 d. It seems 
that a complete, detailed model could be found to explain this 
system along these lines. 

Bot h components in HE 1414-0848 are DA white dwarfs 
( Nap iwotzki et al 120021) . as is the secondary of PG 1115+116. 
The hydrogen in the spectra of these stars suggests that the sur- 
face layer that formed after the envelope was ejected is still 
present. However, the primary i n PG 1115+116 is a DB white 
dwarf. As lMaxted et al.l ( 2002a) point out, the giant phase of a 
helium star could be the explanation for this and the scenario 
sketched above might indeed describe the formation of this sys- 
tem. 

8. Conclusions 

We investigated several formation scenarios for the observed 
ten double white dwarfs listed in Tableland present the best 
models in Table|5] We draw five main conclusions: 

- The scenario where the first mass-transfer phase is stable 
and conservative, followed by a common envelope with 
spiral-in based on energy conservation (see Eq.[4} cannot 
explain the observed masses and periods of all double white 
dwarfs. 

- The scenario with envelope ejection based on angular- 
momentum balance followed by ejection of the second en- 
velope with either energy or angular-momentum balance 
can explain the observed masses and orbital periods very 
well. 

- Including the age difference as a quality criterion for model 
solutions produces strong restrictions to the selection of so- 
lutions and makes it much more difficult to find acceptable 
solutions. 

- By taking into account the possibilities that mass is lost 
either from the donor or from the accretor, we show that 
the formation of the close double white dwarfs can be ex- 
plained if the mass carries the specific angular momentum 
of one of the two binary members. 

- In particular, the scenario in which the envelope of the pri- 
mary is lost with the specific orbital angular momentum of 
the donor star, followed by envelope loss of the secondary 
with the specific orbital angular momentum of its compan- 
ion star (ydTa) can explain the observed masses, orbital pe- 
riods and age differences well. 
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215 


0.0 


1.52 


1.93 


1.83 


1.05 


131.8 


4.69 


43.20 


0.39 


0.34 


0.14 


42 


1101 


7dffcc 


1.10 


1.06 


215 


1059 


392. 


1.48 


2.87 


1.86 


1.54 


21.45 


4.75 


41.47 


0.39 


0.34 


0.14 


43 


1101 


7d7a 


1.11 


0.98 


215 


268 


25. 


0.69 


2.87 


2.40 


1.20 


22.08 


6.13 


24.87 


0.39 


0.34 


0.14 


44 


1101 


7d7a 


1.10 


0.98 


215 


322 


50. 


0.74 


2.87 


2.34 


1.23 


22.02 


5.97 


28.23 


0.39 


0.34 


0.14 


45 


1101 


7d7a 


1.10 


0.95 


215 


363 


69. 


1.10 


2.34 


2.03 


1.15 


5.150 


6.79 


5.880 


0.30 


0.26 


0.14 



Table 6. Selected model solutions for the double envelope-ejection scenario. This table is the total list of which an excerpt 
can be found as Table|5]in the article. The first nine columns show the number of the entry, the double white dwarf that the 
model is a solution for, the prescription used, the two envelope-ejection parameters, the age difference of the two components 
as observed and in the model (At) in Myr, the relative difference between the observed and model age difference, defined as 
A(At) = j AT TOy At ?i» I m %^ the time of the formation of the double white dwarf since the ZAMS (t 2 ) in Gyr. The last nine columns 
list binary parameters: the initial (ZAMS) masses, mass ratio and orbital period, the intermediate mass ratio and period and the 
final masses and period. (The table is continued on the next page) 
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Nr. 


WD 


Presc. 


7i 


72, 


At (Myr) 


A(At) 


T 2 




M 2i 


In 




<?2m 


Pm 




M 2 f 


P f 












Obs 


Mdl 


% 


Gyr 


M e 


M e 




d 




d 


M e 


Me 


d 


46 


1115 


7s«ce 


1.79 


1.00 


160 


239 


49. 


0.50 


3.70 


2.94 


1.26 


1693. 


3.58 


980.4 


0.82 


0.69 


30.09 


47 


1115 


TsQ'cc 


1.95 


1.00 


160 


203 


27. 


0.73 


2.90 


2.59 


1.12 


2088. 


3.24 


1017. 


0.80 


0.67 


30.09 


48 


1115 


7s»ce 


1.90 


1.00 




1UJ 


2.9 


0.50 


3.38 


2.94 


1.15 


1960. 


3.58 


980.4 


0.82 


0.69 


30.09 


49 


1115 


7s7s 


1.80 


1.53 


160 


190 


19. 


0.42 


3.89 


3.13 


1.24 


2284. 


3.51 


1314. 


0.89 


0.75 


30.09 


50 


1115 


TsTs 


1.73 


1.62 


160 


313 


95. 


0.56 


3.74 


2.83 


1.32 


2340. 


3.17 


1500. 


0.89 


0.75 


30.09 


51 


1115 


7a<*ce 


1.41 


0.93 


160 


240 


50. 


0.32 


5.42 


3.42 


1.58 


201.2 


3.84 


1012. 


0.89 


0.75 


30.09 


52 


1115 




1.00 


0.98 


160 


831 


419. 


0.94 


4.89 


2.37 


2.07 


153.1 


3.01 


1053. 


0.79 


0.66 


30.09 


53 


1115 


7a7a 


1.49 


1.03 


160 


229 


43. 


0.31 


5.42 


3.47 


1.56 


201.5 


3.88 


832.5 


0.89 


0.75 


30.09 


54 


1115 


7a7a 


1.08 


1.08 


160 


546 


241. 


0.65 


4.89 


2.69 


1.82 


87.43 


3.47 


814.0 


0.77 


0.65 


30.09 


55 


1115 


7d«ce 


0.97 


0.93 


160 


240 


50. 


0.32 


5.42 


3.42 


1.58 


201.2 


3.84 


1012. 


0.89 


0.75 


30.09 


56 


1115 


7d7a 


1.32 


1.00 


1UU 


164. 


2.5 


0.46 


3.51 


3.02 


1.17 


1378. 


3.90 


427.6 


0.77 


0.65 


30.09 


57 


1115 


7d7a 


0.97 


1.04 


160 


240 


50. 


0.32 


5.42 


3.42 


1.58 


201.2 


3.84 


1012. 


0.89 


0.75 


30.09 


58 


1204 


7s«ce 


2.15 


0.83 


Ou 




70. 


5.37 


1.27 


1.25 


1.01 


630.3 


2.65 


274.3 


0.47 


0.41 


1.60 


59 


1204 


7s*ce 


1.05 


0.96 


80 


225 


181. 


1.58 


2.06 


1.81 


1.14 


32.89 


3.58 


256.2 


0.51 


0.44 


1.60 


60 


1204 


7s7s 


1.89 


1.25 


ou 


ou 


0.5 


0.31 


3.60 


3.21 


1.12 


69.10 


6.07 


39.45 


0.53 


0.46 


1.60 


61 


1204 


'Y'd^ce 


1.45 


0.96 


80 


225 


181. 


1.58 


2.06 


1.81 


1.14 


32.89 


3.58 


256.2 


0.51 


0.44 


1.60 


62 


1204 




1.29 


0.87 


ou 


4.Q8 


522. 


1.85 


2.06 


1.72 


1.20 


32.72 


3.40 


277.6 


0.51 


0.44 


1.60 


63 


1204 


7a7a 


2.42 


0.94 


80 


143 


78. 


0.42 


3.34 


2.87 


1.17 


15.40 


6.08 


37.29 


0.47 


0.41 


1.60 


64 


1204 


7a7a 


1.06 


1.08 


80 


1387 


1634. 


1.66 


3.34 


1.79 


1.87 


14.56 


3.79 


172.9 


0.47 


0.41 


1.60 


65 


1204 


7dC C c 


0.72 


0.96 


ou 


225 


181. 


1.58 


2.06 


1.81 


1.14 


32.89 


3.58 


256.2 


0.51 


0.44 


1.60 


66 


1204 


7d«ce 


1.47 


0.87 


80 


138 


72. 


4.66 


1.32 


1.30 


1.01 


606.6 


2.76 


263.4 


0.47 


0.41 


1.60 


67 


1204 


7d7a 


1.10 


0.95 


80 


74 


7.8 


0.26 


3.89 


3.42 


1.14 


38.82 


5.96 


51.85 


0.57 


0.50 


1.60 


68 


1204 


7d7a 


1.09 


0.92 


80 


100 


25. 


0.38 


3.34 


2.98 


1.12 


15.47 


6.32 


19.99 


0.47 


0.41 


1.60 


69 


1349 




1.56 


1.05 


o 


129 


0.0 


2.17 


1.68 


1.63 


1.03 


108.7 


4.48 


352.9 


0.37 


0.46 


2.21 


70 


1349 




1.45 


1.01 


o 


461 


0.0 


4.52 


1.35 


1.32 


1.03 


105.9 


3.77 


364.5 


0.35 


0.44 


2.21 


71 


1349 


7s7s 


1.56 


1.56 


o 


129 


0.0 


2.17 


1.68 


1.63 


1.03 


93.59 


4.58 


317.3 


0.36 


0.45 


2.21 


72 


1349 


7s7s 


1.52 


1.63 


o 


384 


0.0 


2.55 


1.63 


1.55 


1.05 


112.7 


4.25 


379.4 


0.37 


0.46 


2.21 


73 


1349 




1.01 


1.02 





1776 


0.0 


2.26 


2.72 


1.61 


1.69 


14.85 


4.42 


360.4 


0.37 


0.46 


2.21 


74 


1349 


7a7a 


1.35 


0.97 





975 


0.0 


1.53 


2.59 


1.83 


1.41 


10.16 


5.38 


206.9 


0.34 


0.43 


2.21 


75 


1349 


7dQ' C c 


0.97 


1.05 


o 


129 


0.0 


2.17 


1.68 


1.63 


1.03 


108.7 


4.48 


352.9 


0.37 


0.46 


2.21 


76 


1349 


7d«ce 


0.96 


1.02 


o 


215 


0.0 


2.26 


1.68 


1.61 


1.04 


108.5 


4.42 


360.4 


0.37 


0.46 


2.21 


77 


1349 


7d7a 


0.95 


0.98 


o 


101 


0.0 


1.58 


1.86 


1.81 


1.03 


63.44 


5.19 


241.2 


0.35 


0.44 


2.21 


78 


1349 


7d7a 


0.97 


1.03 


o 


235 


0.0 


2.17 


1.70 


1.63 


1.04 


106.4 


4.48 


352.9 


0.37 


0.46 


2.21 


79 


1414 


7s»ce 


1.52 


0.71 


200 


188 


5.9 


0.43 


3.51 


3.09 


1.14 


70.81 


5.99 


358.3 


0.52 


0.66 


0.52 


80 


1414 


7s7s 


1.49 


1.76 


200 


112 


44. 


0.87 


2.55 


2.43 


1.05 


525.4 


4.09 


1783. 


0.59 


0.76 


0.52 


81 


1414 




2.13 


0.71 


200 


196 


1.9 


0.43 


3.56 


3.09 


1.15 


51.80 


5.99 


358.3 


0.52 


0.66 


0.52 


82 


1414 


7a7a 


1.10 


1.04 


200 


290 


45. 


0.47 


3.99 


3.02 


1.32 


41.73 


5.08 


1363. 


0.59 


0.76 


0.52 


83 


1414 


7a7a 


1.04 


1.04 


200 


325 


63. 


0.50 


3.99 


2.94 


1.36 


41.61 


4.95 


1407. 


0.59 


0.76 


0.52 


84 


1414 


7d*ce 


0.95 


0.71 


200 


188 


5.9 


0.43 


3.51 


3.09 


1.14 


70.81 


5.99 


358.3 


0.52 


0.66 


0.52 


85 


1414 


7d7a 


0.95 


0.99 


200 


188 


5.9 


0.43 


3.51 


3.09 


1.14 


70.81 


5.99 


358.3 


0.52 


0.66 


0.52 


86 


1414 


7d7a 


0.90 


1.00 


200 


167 


17. 


0.39 


3.65 


3.21 


1.14 


75.14 


5.96 


511.9 


0.54 


0.69 


0.52 


87 


1704a 


7s fee 


1.67 


0.60 


-20 


52 


360. 


1.41 


2.06 


1.88 


1.09 


40.37 


3.66 


65.66 


0.51 


0.36 


0.14 


88 


1704a 


7s fee 


1.88 


0.62 


-20 


15 


175. 


1.17 


2.19 


2.00 


1.09 


93.52 


3.79 


66.89 


0.53 


0.37 


0.14 


89 


1704a 


7s fee 


2.05 


0.43 


-20 


7 


135. 


1.36 


2.03 


1.90 


1.07 


252.8 


3.51 


96.02 


0.54 


0.38 


0.14 


90 


1704a 


7s7s 


1.67 


1.52 


z,u 




360. 


1.41 


2.06 


1.88 


1.09 


40.37 


3.66 


65.66 


0.51 


0.36 


0.14 


91 


1704a 


7s7s 


1.88 


1.50 


-20 


15 


175. 


1.17 


2.19 


2.00 


1.09 


93.52 


3.79 


66.89 


0.53 


0.37 


0.14 


92 


1704a 


7aCce 


2.64 


0.59 


-20 


23 


215. 


1.22 


2.16 


1.98 


1.09 


87.04 


3.74 


69.77 


0.53 


0.37 


0.14 


93 


1704a 


7a«ce 


1.90 


0.43 


z,u 


4.QQ 


2596. 


1.36 


2.43 


1.90 


1.27 


43.81 


3.51 


96.02 


0.54 


0.38 


0.14 


94 


1704a 


7a7a 


1.57 


1.09 


-20 


884 


4521. 


1.12 


3.56 


2.03 


1.75 


43.52 


3.95 


50.36 


0.51 


0.36 


0.14 


95 


1704a 


7dQ'ce 


1.11 


0.60 


-20 


52 


360. 


1.41 


2.06 


1.88 


1.09 


40.37 


3.66 


65.66 


0.51 


0.36 


0.14 


96 


1704a 


7dQ'ce 


1.24 


0.62 


-20 


15 


175. 


1.17 


2.19 


2.00 


1.09 


93.52 


3.79 


66.89 


0.53 


0.37 


0.14 


97 


1704a 


7dQ'ce 


1.37 


0.34 


-20 


2 


110. 


1.48 


1.93 


1.86 


1.04 


294.5 


3.33 


120.2 


0.56 


0.39 


0.14 


98 


1704a 


7d7a 


1.11 


1.13 


-20 


52 


360. 


1.41 


2.06 


1.88 


1.09 


40.37 


3.66 


65.66 


0.51 


0.36 


0.14 


99 


1704a 


7d7a 


1.24 


1.11 


-20 


15 


175. 


1.17 


2.19 


2.00 


1.09 


93.52 


3.79 


66.89 


0.53 


0.37 


0.14 
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Nr. 


WD 


Presc. 


7i 


72, 


At (Myr) 


A(At) 


1"2 


M n 


M 2i 




Pi 


<?2m 


Pm 




M 2f 


Pt 










Q'ce2 


Obs 


Mdl 


% 


Gyr 


M e 


M 




d 




d 


Me 


Me 


d 


100 


1704b 


7s fee 


1.65 


0.53 


20 


292 


1360. 


0.73 


2.83 


2.59 


1.09 


47.21 


6.37 


161.8 


0.41 


0.58 


0.14 


101 


1704b 




1.74 


0.76 


20 


285 


1326. 


0.75 


2.76 


2.55 


1.08 


49.12 


6.40 


107.3 


0.40 


0.57 


0.14 


102 


1704b 


7s7s 


1.63 


1.43 


Z.U 




927. 


0.58 


2.98 


2.79 


1.07 


50.44 


6.54 


208.9 


0.43 


0.61 


0.14 


103 


1704b 


7a«ce 


2.36 


0.46 


20 


231 


1054. 


0.58 


3.05 


2.79 


1.09 


33.85 


6.54 


208.9 


0.43 


0.61 


0.14 


104 


1704b 


7aCce 


1.59 


0.47 


20 


1017 


4985. 


1.45 


2.83 


1.95 


1.45 


20.80 


5.08 


158.7 


0.38 


0.55 


0.14 


105 


1704b 


7a7a 


1.63 


1.01 


20 


498 


2391. 


0.86 


3.02 


2.43 


1.24 


31.86 


5.78 


527.9 


0.42 


0.60 


0.14 


106 


1704b 


7a7a 


2.56 


0.99 


20 


205 


927. 


0.58 


2.98 


2.79 


1.07 


50.44 


6.54 


208.9 


0.43 


0.61 


0.14 


107 


1704b 


7d Q 'cc 


1.03 


0.15 


Z.U 


182 


810. 


2.23 


1.68 


1.65 


1.01 


212.1 


4.08 


478.6 


0.41 


0.58 


0.14 


108 


1704b 


7d«ce 


1.00 


0.76 


20 


332 


1562. 


0.75 


2.87 


2.55 


1.12 


33.29 


6.40 


107.3 


0.40 


0.57 


0.14 


109 


1704b 


7d7a 


0.97 


0.99 


20 


205 


927. 


0.58 


2.98 


2.79 


1.07 


50.44 


6.54 


208.9 


0.43 


0.61 


0.14 


110 


1704b 


7d7a 


0.99 


0.99 


Z.U 


292 


1360. 


0.73 


2.83 


2.59 


1.09 


47.21 


6.37 


161.8 


0.41 


0.58 


0.14 


111 


2209 


7s<2ce 


1.69 


0.54 


^ on 


^17 

J L l 


3.3 


1.45 


2.37 


1.95 


1.21 


148.5 


3.55 


168.2 


0.55 


0.55 


0.28 


112 


2209 


7s«ce 


1.56 


0.88 


JUU 




10. 


0.75 


3.79 


2.55 


1.49 


87.01 


4.48 


113.8 


0.57 


0.57 


0.28 


113 


2209 


7s7s 


1.62 


1.69 


500 


477 


4.5 


1.37 


2.40 


2.00 


1.20 


165.4 


3.64 


258.7 


0.55 


0.55 


0.28 


114 


2209 


7s7s 


1.62 


1.63 


JUU 


262 


48. 


1.20 


2.37 


2.16 


1.09 


150.0 


3.93 


304.6 


0.55 


0.55 


0.28 


115 


2209 


7a«ce 


1.67 


0.70 


500 


690 


38. 


0.90 


3.74 


2.40 


1.56 


58.62 


4.36 


138.0 


0.55 


0.55 


0.28 


116 


2209 


7a«ce 


1.30 


0.54 


5fin 

JUU 


1245 


149. 


1.45 


3.74 


1.95 


1.92 


57.09 


3.55 


168.2 


0.55 


0.55 


0.28 


117 


2209 


7a7a 


1.38 


1.10 


500 


542 


8.4 


0.70 


4.15 


2.62 


1.58 


97.10 


4.16 


677.7 


0.63 


0.63 


0.28 


118 


2209 


7a7a 


1.14 


1.16 


500 


846 


69. 


1.00 


4.15 


2.31 


1.80 


95.54 


3.66 


821.2 


0.63 


0.63 


0.28 


119 


2209 


7dC C c 


1.06 


0.53 


500 


347 


31. 


1.35 


2.31 


2.03 


1.14 


80.61 


3.76 


169.5 


0.54 


0.54 


0.28 


120 


2209 


7dCcc 


1.12 


0.88 


500 


559 


12. 


0.75 


3.84 


2.55 


1.50 


71.51 


4.48 


113.8 


0.57 


0.57 


0.28 


121 


2209 


7d7a 


1.06 


1.05 


500 


494 


1.3 


0.67 


3.94 


2.65 


1.48 


83.11 


4.50 


229.2 


0.59 


0.59 


0.28 


122 


2209 


7d7a 


1.04 


1.05 


500 


340 


32. 


0.50 


4.15 


2.94 


1.41 


98.45 


4.67 


294.3 


0.63 


0.63 


0.28 
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